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PROTEIN ENGINEERING OF GLUCO AMYLASE TO INCREASE 
PH OPTIMUM, SUBSTRATE SPECIFICITY AN D THERMOSTABILITY 

FIELD OF THE INVENTION 

5 The field of the invention relates to mutations to 

produce a fungal glucoamylase enzyme that is more 
selective for the production of glucose rather than the 
a-1,6 linked disaccharide isomaltose, is more 
thermostable, and has increased pH optimum and produces 
10 increased amounts of glucose compared to wildtype 
enzymes . 

BACKGROUND OF THE INVENTION 

Glucoamylase (EC 3.2.1.3) is a carbohydrase . Dis- 

15 covered in 1951, it is an exo-hydrolase that cleaves d- 
glucose from the nonreducing ends of 
maltooligosaccharides, attacking a- (1,4)-, and at a 
much slower rate, or- (1, 6) -glucosidic bonds. It is one 
of more than one hundred carbohydrases (EC 3.2.1) that 

20 cleave O-glycosidic bonds of either a- or 0- 

conf iguration. The functional and structural related- 
ness of these enzymes is reflected in the presence of 
at least three discrete regions of sequence homology 
between glucoamylase and several a-amylases, <*- 

25 glucosidases, and transglucanosylases [Svensson, 1988], 
and a similar domain structure to carbohydrases that 
attack insoluble substrates [Knowles et al., 1987; 
Svensson et al., 1989)]. Aspergillus awamori 
glucoamylase (1, 4-a-D-glucan glucohydrolase; EC 

30 3.2.1.3) 'is one of the most important of the 
glucoamylases. 

Glucoamylase is primarily used in industry for the 
production of high-f ructose corn sweeteners in a proc- 
ess that involves 1) a-amylase to hydrolyze starch to 

35 maltooligosaccharides of moderate length (dextrin) ; 2) 
Glucoamylase to hydrolyze dextrin to glucose; and 3) 
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glucose isomerase to convert glucose to fructose . Corn 
sweeteners have captured over 50% of the U. S. 
sweetener market, and the three enzymes used to make 
them are among the enzymes made in highest volume. In 
5 addition, glucose produced by glucoamylase can be 

crystallized or used in fermentation to produce organic 
products such as citric acid, ascorbic acid, lysine, 
glutamic acid or ethanol for beverages and fuel. 
Approximately 12% of the country's corn production is 

10 processed with glucoamylase. Although glucoamylase has 
been successfully used for many years, it would be a 
more attractive product if it produced higher amounts 
of glucose instead of disaccharides, if it were more 
stable, and if it could be used in the same vessel with 

15 glucose isomerase* 

Glucoamylase does not give 100% yield of glucose 
from dextrin because it makes various di- and 
trisaccharides, especially isomaltose and 
isomaltotriose, from glucose [Nikolov et al., 1989]. 

20 These products, formed at high substrate 

concentrations, result from the ability of glucoamylase 
to form a- (1, 6) -glucosidic bonds. Glucoamylase is not 
as thermostable as either or- amylase or glucose 
isomerase. The optimum pH of GA (pH4-4.5) is lower 

25 than that of aamylase (pH5.5-6.5) and glucose isomerase 
(pH7-8) . Therefore glucoamylase hydrolysis must be 
done separately from the other enzymatic reactions in a 
different vessel and at lower temperatures, causing 
higher capital costs. 

30 Glucoamylase from the filamentous fungus 

Aspergillus niger is the most widely used glucoamylase, 
and its biochemical properties have been extensively 
characterized. This enzyme is found mainly in two 
forms, GAI (616 amino acids; referred to as AA 

35 hereinafter) and GAII (512 AA) , differing by the pres- 
ence in GAI of a 104-AA C-terrninal domain required for 
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adsorption to native starch granules [Svensson et al., 
1982; Svensson et al., 1989]. Both forms have a 
catalytic domain (AA1-440) followed by a Ser/Thr-rich, 
highly O-glycosylated region (AA44 1-512) [Gunnarsson et 
5 al., 1984]. The first thirty residues of this region 
are included in the three-dimensional structure of the 
enzyme [Aleshin et al., 1994; 1996; Stoffer et al., 
1995] ; they wrap around the catalytic domain like a 
belt. There is strong AA sequence homology among 

10 fungal glucoamylase' s in four distinct regions of the 

catalytic domain that correspond to the loops that form 
the substrate binding site [Itoh et al., 1987] . In A. 
niger glucoamylase these regions are AA35-59, AA104- 
134, AA162-196, and AA300-320. The second and third 

15 regions partially or completely overlap the three 

regions of homology to ct -amylases [Svensson, 1988] . In 
addition, the raw starch binding domain (AA512-616) has 
high homology to similar domains from several starch- 
degrading enzymes [Svensson et al. f 1989]. 

20 Kinetic analysis showed that the substrate binding 

site is composed of up to seven subsites [Savel'ev et 
al., 1982] with hydrolysis occurring between subsites 1 
and 2. The pK a ' s of hydrolysis, 2.75 and 5.55 [Savel'ev 
and Firsov, 1982] , suggest that carboxylic acid 

25 residues at subsites 1 and 2 provide the catalytic acid 
and base for hydrolysis. Chemical modification 
experiments showed that three highly conserved 
residues, Aspl76, Glul79, and Glul80, are protected and 
are in the active site, suggesting that one or more of 

30 them are the possible catalytic residues [Svensson et 

al., 1990] . Chemical modification experiments also in- 
dicated that the highly conserved residue Trpl20 is 
essential, and is located in subsite 4 [Clarke and 
Svensson, 1984] . Trpl20 is homologous to Trp83 of 

35 Aspergillus oyzae a-amylase [Clarke and Svensson, 

1984] , which is also located in the active site of that 
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enzyme [Matsuura et al., 1984] . Site directed 
mutagenesis studies have indicated that Glul79 is the 
catalytic acid residue, while Glu400 is the catalytic 
base residue [Frandse et al, 1994; Harris et al, 1993; 
5 Sierks et al, 1990] 

Glucoamylases from A. niger [Svensson et al., 
1983; Boel et a J . , 1984] and Aspergillus awamori 
[Nunberg et al., 1984] have been cloned and sequenced, 
and have identical primary structures. Innis et al- 
io [1985] and more recently Cole et al. [1988] have 

developed vectors (pGAC9 and pPM18, respectively) for 
glucoamylase expression in yeast, allowing convenient 
manipulation and testing of glucoamylase mutants. 

15 SUMMARY OF THE INVENTION 

According to the present invention, a fungal 
glucoamylase (1,4-a-D-glucan glucohydrolase ; EC 3.2.1) 
with decreased thermal inactivation (increased 

20 thermostability) and reduced isomaltose formation 

provided by the mutation Asn20Cys coupled with Ala27Cys 
forming a disulfide bond between the two is provided. 
Cumulative thermostability is also provided for GA by 
including the mutation Asn20Cys coupled with Ala27Cys 

25 and at least one mutation from Table 13. An 

engineered GA including Ser30Pro, Glyl37Ala, and 
Asn20Cys coupled with Ala27Cys provides even more 
thermostability. Cumulative thermostability is also 
provided for GA by including the mutation Asn20Cys 

30 coupled with Ala27Cys and at least two mutations from 
Table 13 . 

The present invention also provides a fungal 
glucoamylase with reduced isomaltose formation 
including an Asn20Cys coupled with Ala27Cys mutation 
35 (S~S mutation) and at least one mutation selected from 
Table 14. In an embodiment Asn20Cys coupled with 
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Ala27Cys mutation and a 311-314Loop (also referred to 
as 300Loop) mutation are included in an engineered GA. 
In a further preferred embodiment the engineered 
glucoamylase with reduced isomaltose formation includes 
5 Asn20Cys coupled with Ala27Cys mutations Ser30Pro and 
Glyl37Ala. 

The present invention also provides engineered 
fungal glucoamylase including a 311-3l4Loop mutation 
whereby reduced isomaltose formation is provided by the 

10 mutation. In a further embodiment fungal glucoamylase 
including a 311-314Loop mutation and at least one 
mutation from Table 14 are prepared whereby cumulative 
reduced isomaltose formation is provided by the 
additional mutation. 

15 The present invention provides a fungal 

glucoamylase including a mutation Ser411Ala whereby 
increased pH optimum and reduced isomaltose formation 
is provided by the mutation. In an embodiment the 
Ser411Ala mutation is combined with at least one 

20 mutation from Table 15 whereby cumulative increased pH 
optimum is provided by the mutations. In an embodiment 
the Ser411Ala mutation is combined with at least one 
mutation from Table 14 whereby cumulative reduced 
isomaltose formation is provided by the mutations. 

25 In a further embodiment an engineered fungal 

glucoamylase includes a mutation Ser411Ala and a 
mutation pair Asn20Cys coupled with Ala27Cys forming a 
disulfide bond between the two members of the pair 
whereby increased thermal stability, increased pH 

30 optimum and reduced isomaltose formation are provided 
by the mutations. 

In a still further embodiment a fungal 
glucoamylase is engineered to include a Ser411Ala 
mutation and a mutation pair Asn20Cys coupled with 

35 Ala27Cys forming a disulfide bond between the two 

members of the pair and a 3 11 -3 14 Loop mutation whereby 
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increased thermal stability, increased pH optimum and 
reduced isomaltose formation are provided by the 
mutations. 

The present invention provides a method to obtain 
5 a fungal glucoamylase with reduced isomaltose formation 
by designing mutations to decrease the a- (1,6)- . 
glucosidic linkage affinity of GA. 

The present invention also provides a method to 
obtain a fungal glucoamylase with decreased thermal 
10 inactivation by designing mutations to decrease the 
enzyme's conformational entropy of unfolding and/or 
increase stability of or-helices, increase disulfide 
bonds, hydrogen bonding, electrostatic interactions, 
hydrophic interactions, Vanderwalls interactions and 
15 packing compactness. 

The present invention also provides a fungal 
glucoamylase with increased pH optimum including 
changing the polarity, charge distribution and hydrogen 
bonding in the microenvironment of the catalytic base 
20 Glu400. 

The present invention also provides a method of 
genetically engineering glucoamylase carrying at least 
two cumulatively additive mutations. Individual 
mutations are generated by site-directed mutagenesis. 

25 These individual mutations are screened and those 

selected which show increased pH optimum and which show 
decreased irreversible thermal inactivation rates or 
reduced isomaltose formation. Site directed 
mutagenesis is then performed to produce enzymes 

30 carrying at least two of the isolated selected 

mutations. Finally the engineered enzymes are screened 
for cumulatively additive effects of the mutations on 
thermal stabilizing or reduced isomaltose formation by 
the produced enzymes carrying at least two of the 

35 isolated selected mutations. Alternatively, the 

engineered enzyme is screened for cumulatively additive 



WO 98/03639 



PCT/US97/12983 



effects of both of the mutations on pH optimum, 
thermostability and/or reduced isomaltose formation by 
the produced enzymes carrying at least two of the 
isolated selected mutations. 
5 Vectors for each of the mutations and mutation 

combinations are also provided by the present invention 
as well as host cells transformed by the vectors . 

DESCRIPTION OF THE DRAWINGS 

10 

Other advantages of the present invention will be 
readily appreciated as the same becomes better 
understood by reference to the following detailed 
description when considered in connection with the 

15 accompanying drawings wherein: 

FIGURE 1 is a graph showing the relationship 
between temperature and kd for wild- type (•) and 
proline substituted mutant GA's: S30P (■) , D345P (v) , 
E408P (O) in Example 1. 

20 FIGURE 2 is a graph showing effect of temperature 

on first -order thermoinactivation rate coefficients of 
wild-type (O) , A27C (•) , N20C (v) , A27C/N20C (t) , 
A471C/T72C {□) , A27C/N20C/G137A (■) , A27C/N20C/S436P 
(O) AND G137A/S436P (♦) glucoamylases measured in pH 

25 4.5 buffer. 

FIGURE 3 is a graph showing initial reaction rates 
of wild-type (O) # A27C/N20C (•) , A471C/T72C (v) and 
A29C/N20C/G137A (▼) glucoamylases with 4% maltose in 
0.05 M sodium acetate (pH 4.5) as substrate at 

30 temperatures from 60°C to 76°C. 

FIGURE 4 is a graph showing the effect of 
temperature on the activity of wildtype and mutant GA. 
Error bars represent the standard deviation from three 
assays. Wildtype (•) , S30P/G137A (□) , S-S/S30P/G137A 

35 (a). 

FIGURE 5A-C are graphs showing the effect of 
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temperature on irreversible thermal inactivation rate 
coefficients of wildtype and mutant GA. Fig. 5A 
Wildtype (•) , S30P (■) , G137A (a), S30P/G137A (□) ; Fig. 
5B Wildtype (•) , S30P (■) , S-S (hexagon) , S-S/S30P 
5 (filed circle with empty center); Fig. 5C Wildtype (•) , 
S30P/G137A (□) , S-S/S30P (filed circle with empty 
center) , S-S/S30P/G137A (a) . 

FIGURE 6A-B are graphs showing saccharif ication of 
28% (w/v) Maltrin M100 by wildtype (•) , S30P/G137A (□) 

10 and S-S/S30P/G137A (a) . 

FIGURE 7 is a graph showing the 30% DE 10 
maltodextrin saccharif ication of wildtype (♦) and 
mutant glycoamylases : 300Loop (■) , S30P/G137A (a), S-S 
(•) , S30P/G137A/300Loop (x) , S-S/300Loop (a) , at 55°C f 

15 enzyme concentration was 166.67 /ig/mL in each reaction. 
FIGURE 8 is a graph showing production of 
isomaltose by wildtype (•) and mutant glucoamylases : 
Y116W <■), Y175F (a), R241K (▼) , S411A (♦) , S411G 
(hexagon) , during glucose condensation at 55°C with 30% 

20 (w/v) D-glucose in 0.05M sodium acetate buffer at pH4.4 
with 0.02% sodium azide for 12 days. 

FIGURE 9 is a graph showing the production of 
glucose by wildtype (•) and mutant glucoamylases: Y116W 
(■), Y175F (a), R241K (▼) , S411A (♦) , S411G (hexagon), 

25 during hydrolysis of DE 10 maltodextrin at 55°C with 
28% (w/v) maltodextrin in 0.05M sodium acetate buffer 
at pH4.4 with 0.02% sodium azide for 12 days. 

FIGURE 10 is a graph showing the initial rates of 
glucose production by wildtype (•) and S411A (■) 

30 glucoamylases during DE 10 maltodextrin hydrolysis at 
different pH values. Hydrolysis was performed at 36°C 
with 28% (w/v) maltodextrin in 25mM citrate -phosphate 
buffer at indicated pHs with 0.02% sodium azide for 4 
days. 



35 



WO 98/03639 



PCT/US97/12983 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

The present invention provides mutations for 
increased thermal stability, increased pH optimum and 
5 reduced isomaltose formation in the glucoamylase from 
fungal species which may provide increased glucose 
yields compared to wildtype glucoamylase. Predicted 
structure and known sequences of glucoamylase are 
conserved among the fungal species [Coutino et al, 

10 1994] . As an exemplar Aspergillus awamori glucoamylase 
(1,4-or-D-glucan glucohydrolase ; EC 3.2.1.3; referred to 
as GA herein; SEQ ID No:l) is used, but any other 
fungal species including Aspergillus species 
glucoamylase can be used. The numbering of the 

15 glucoamylase amino acids herein is based on the 
sequence of the exemplar Aspergillus awamori. 
Equivalent amino acid residue numbers are determined 
differently for different fungal species as is known in 
the art [Coutino et al., 1994]. 

20 The present invention provides a fungal 

glucoamylase with decreased thermal inactivation 
(increased thermostability) and decreased isomaltose 
formation provided by engineering the inclusion of a 
mutation pair Asn20Cys coupled with Ala27Cys which 

25 forms a disulfide bond between them (this mutation is 
abbreviated as Asn20Cys/Ala27Cys or S-S) . Additional 
mutations providing decreased thermal inactivation are 
set forth in Summary Table 13. 

Cumulative thermostability is also provided for GA 

30 by including at least two of the mutations in the 

enzyme as for example including mutations Ser30Pro and 
Glyl37Ala. Another example is to engineer S-S with 
Asn20Cys/Ala27Cys in the enzyme or to pair Glyl37Ala 
with S-S. Further, combinations of the individual 

35 mutations set forth in Table 13, particularly with S-S 
coupled with Ser30Pro also provide cumulative 



WO 98/03639 



PCT/US97/12983 



-10- 

thermostability. In general two mutation combinations 
are made but triple mutations can also be constructed. 
As for example, an engineered GA including the three 
mutations: Ser30Pro, Glyl37Ala, and Asn20Cys/Ala27Cys 
5 provides even more thermostability. 

By Asn20Cys coupled with Ala27Cys is meant a pair 
of mutations which is abbreviated as n S-S w or 
Asn20Cys/Ala27Cys and between which is formed a 
disulfide bond as described herein in the Examples. In 

10 general, this is referred to as a single mutation since 
both are required to form the disulfide bond. 

By cumulative is generally meant the additive (or 
nearly additive) effects of two or more mutations on 
the parameter of enzyme activity being measured. 

15 The present invention also provides a fungal 

glucoamylase with reduced isomaltose formation and 
increased glucose yield including the Asn20Cys/Ala27Cys 
mutation (S-S mutation) and at least one mutation 
selected from Table 14. In an embodiment the 

20 Asn20Cys/Ala27Cys mutation and the 311-314Loop 

(300Loop) are included in GA. In a further preferred 
embodiment the engineered glucoamylase with reduced 
isomaltose formation includes Asn20Cys/Ala27Cys and 
with mutations Ser30Pro and Glyl37Ala. 

25 In an embodiment a glucoamylase with the 311-114 

loop mutation is constructed to provide reduced 
isomaltose formation. By the 311-314Loop mutation is 
meant an insertional GA mutant with the sequence 
Tyr311-Tyr312-Asn313-Gly314^Tyr311-Asn-Gly-Asn-Gly-Asn- 

30 Ser-Gln-Gly314 (311-314 Loop; SEQ ID No: 2) . 

The present invention provides a fungal 
glucoamylase including a Ser411Ala mutation whereby 
increased pH optimum and reduced isomaltose formation 
is provided by the mutation. In an embodiment the 

35 Ser411Ala mutation is combined with at least one 

mutation from Table 15 whereby cumulative increased pH 
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optimum is provided by the combined mutations. In a 
further embodiment the Ser411Ala mutation is combined 
with at least one mutation from Table 14 whereby- 
cumulative reduced isomaltose formation is provided by 
5 the mutations. 

In a further embodiment an engineered fungal 
glucoamylase includes a Ser411Ala mutation and the 
mutation pair Asn20Cys/Ala27Cys forming a disulfide 
bond between them whereby increased thermal stability, 

10 increased pH optimum and reduced isomaltose formation 
are provided by the mutations. 

In a still further embodiment a fungal 
glucoamylase including a Ser411Ala mutation and a 
mutation pair Asn20Cys coupled with Ala27Cys forming a 

15 disulfide bond between the two members of the pair and 
a 311-314Loop mutation whereby increased thermal 
stability, increased pH optimum and reduced isomaltose 
formation are provided by the combination of mutations. 
Mutations are indicated by the amino acid being 

20 replaced followed by the residue number followed by the 
replacing amino acid. Amino acids are abbreviated 
either with the three letter code or single letter 
code. Mutations are generated using site directed 
mutagenesis as is known in the art. The sequence and 

25 residue number are from the Wildtype (WT) or nonmutant 
enzyme. Biochemical characterization is performed as 
described herein below and in the Examples. The 
Examples provide exemplars of the analysis for an 
individual mutation to determine it's characteristics 

30 and provide exemplars of analysis for combinations of 
mutations to determine if the combination provides a 
cumulative effects . 

By increased thermostability (or decreased thermal 
inactivation) is meant that at temperatures between 

35 65 °C and 77.5°C the mutants are irreversibly 

inactivated at a decreased rate compared to wildtype. 
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The present invention provides a method to obtain 
fungal glucoaraylases with decreased thermal 
inactivation by designing mutations to decrease the 
rate of irreversible thermal inactivation at 
5 temperatures between 65°C and 77.5°C compared to 
wildtype. This is accomplished by designing 
glucoamylases with decreased thermal inactivation by 
designing mutations to decrease the enzyme's 
conformational entropy of unfolding and/or increase 

10 stability of a-helices, increase disulfide bonds, 

hydrogen bonding, electrostatic interactions, hydrophic 
interactions, Vanderwalls interactions and packing 
compactness. 

Basic mechanisms underlying protein 

15 thermostability and factors influencing reversible and 
irreversible thermal inactivation have been studied 
extensively [Argos et al . , 1979; Klibanov, 1983; 
Wasserman, 1984; Ahem and Klibanov, 1985] . Factors 
involved in stabilizing proteins at high temperatures 

20 include 1) disulfide bonds, 2) noncovalent bonds such 
as salt bridges, hydrogen bonding, and hydrophobic 
interactions, and 3) conformational rigidity [Nosoh and 
Sekiguchi, 1988] . The causes of irreversible 
inactivation at high temperatures include 1) 

25 aggregation, 2) the formation of incorrect structures, 
3) the destruction of disulfide bonds, 4) deamidation 
(especially of Asn at Asn-Gly sequences) , and 5) 
cleavage of Asp-X peptide linkages. It is apparent 
that replacement of even one residue can make a large 

3 0 difference in protein thermostability [Matsumura and 

Aiba, 1985] , due to the small increases in free energy 
(20-30 kJ/mol) usually required to stabilize protein 
tertiary structures [Nosoh and Sekiguchi, 1988] . 
Genetic engineering to increase thermostability (or to 

35 decrease irreversible thermoinactivation) of enzymes 

has been successful in several cases [Perry and Wetzel, 
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1984; Imariaka et al. f 1986; Ahearn et al. f 1987]. 
However, the mechanisms that govern thermostability are 
not fully understood, so that amino acid (AA) 
replacements that promote thermostability are not 
5 accurately predicted [Leatherbarrow and Fersht, 1986; 
Nosoh and Sekiguchi, 1988; Pakula and Sauer, 1989]. 
The method of the present invention allows for more 
accurate prediction. 

By increased pH optimum is meant that the enzyme 

10 is functional at a higher pH, above that of wildtype. 

The present invention also provides a method to 
design a fungal glucoamylase with increased pH optimum 
by changing the polarity, charge distribution and 
hydrogen bonding in the microenvironment of the 

15 catalytic base Glu400. For example, mutants S411G and 
S411A were designed to remove the hydrogen bond between 
Ser411 and Glu400 (see Example 8) . 

By increased selectivity is meant that there is 
decreased isomaltose formation due to decrease in the 

20 production of undesirable a- (1, 6) -linked byproducts 
(reversion products) at high glucose concentrations 
[Lee et al . , 1976] . As described above, GA hydrolyzes 
and synthesizes both o?-(l,4) and a- (1,6) glucosidic 
bonds. Increasing selectivity indicates that the 

25 enzyme synthesizes <xl,6 linked products at a lower rate 
than wildtype as shown by reduced levels of isomaltose 
formation in condensation reactions with 30% glucose as 
a substrate compared to wildtype GA. Additionally, 
improved selectivity may result in increased glucose 

30 yields in saccharif ication reactions using 28% DE 10 
maltodextran as a substrate. 

The present invention provides a method to obtain 
a fungal glucoamylase with reduced isomaltose formation 
by designing mutations to decrease the a- (1,6)- 

35 glucosidic linkage affinity of GA. That is mutations 
are designed in the active site to reduce isomaltose 
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formation due to glucose condensation- The mutations 
are designed to have decreased ability to synthesize 
isomaltose while maintaining at least partial wildtype 
ability to digest al-4 linked substrates resulting in a 
lower ratio of isomaltose formation to glucose 
formation than wildtype. These mutations are made at 
positions that are not completely conserved based on 
homology analysis. 

Kinetic studies have indicated that there are five 
to seven glucosyl binding subsites, and the catalytic 
site is located between subsites 1 and 2 [Hiromi et 
al. f 1973, Hiromi et al., 1983, Meagher et al., 1989, 
Tanaka et al., 1983], The solved three-dimensional 
structure of the catalytic domain of glucoamylase from 
Aspergillus awamori var X100, which has about 95% 
homology with the corresponding regions of GAs from 
Aspergillus awamori and Aspergillus niger [Coutinho & 
Reilly, 1994] , contains thirteen alpha-helices, twelve 
of which are arranged in pairs forming an alpha/alpha 
barrel [Aleshin et al., 1992, Aleshin et al., 1994]. 
The active site is located in the cavity of the barrel 
center. In addition, homology analysis of thirteen 
amino acid sequences of glucoamylases showed that five 
conserved regions define the active site [Coutinho & 
Reilly, 1994] . The mechanism of GA catalysis involves 
two carboxyl groups [Hiromi et al., 1966], Glul79 and 
Glu400 (in Aspergillus awamori or Aspergillus niger) 
[Frandsen et al., 1994, Harris et al., 1993, Sierks et 
al., 1990]. Glul79 protonates the oxygen in the 
glycosidic linkage, acting as general acid acatalyst, 
and Glu400 activates water (WatBOO) for nucleophilic 
attack at carbon Ol, acting as a general base catalyst 
[Frandsen et al., 1994]. The crystal structures of 
glucoamylase complexed with the pseudotetrasaccharides 
(acarbose and D-gluco-dihydroacarbose) , showed that 
there are two different binding conformers, pH 4 -type 
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and pH 6 -type, for pseudotetrasaccharides at pH 4 
[Aleshin et al. # 1996, Stoffer et al. f 1995]. Binding 
of the first two sugar residues of the 
pseudotetrasaccharides is the same, but there is an 
5 extraordinary variation in binding of the third and 
fourth sugar residues of the pseudotetrasaccharides 
[Stoffer et al., 1995] . 

The substrate specificity of an enzyme is 
determined by its ability to form a stable complex with 

10 a ligand in both the ground state and the transition 
state. The stability of the enzyme-ligand complex is 
affected by steric constraints, hydrogen bonding, van 
der Waal's and electrostatic forces, and hydrophobic 
contacts [see generally Fersht, 1985 Enzyme Structure 

15 and Mechanism, 2 nd edition, Freeman, San Francisco] . 
Site-directed mutagenesis was used to construct 
mutations at residues 119 and 121 to alter the hydrogen 
bonding between enzyme and substrate. Atom OG of 
Serll9 hydrogen bonds to the 3 -OH of the fourth sugar 

20 residue of pseudo- tetrasaccharides only in the pH 6- 
type conformer, whereas the amide nitrogen of Glyl21 
hydrogen bonds to the 6 -OH of the third sugar residue 
in both pH 4 -type and pH 6 -type conf ormers . These 
mutations are designed to change substrate specificity 

25 (decrease alpha- 1,6 condensation reactions) while 

maintaining wild- type ability to hydrolyze alpha-1,4 
linked substrates. Serll9 is not conserved and is 
replaced by Ala, Pro and Glu in other GAs . Mutant 
S119E was designed to strengthen the hydrogen bond 

30 between the enzyme and the fourth sugar residue of the 
substrate to stabilize the pH 6 -type conformer, and to 
bring a negative charge near subsite 4 in order to 
increase electrostatic interactions in active site. 
Mutant S119G was designed to remove the same hydrogen 

3 5 bond in order to destabilize the pH 6 -type conformer. 
Mutant S119W was designed to remove the same hydrogen 
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bond and to increase the hydrophobic interactions 
between the enzyme and the pH 6 -type conformer. Glyl21 
is highly conserved in all glucoamylase sequences 
except in Clostridium sp. G005 GA, which has high a-1,6 
5 activity and in which Gly is replaced by Thr. Since 
the <f> and ^ angles of Glyl21 would allow an alanine in 
this position without causing a conformation 
distortion, G121A was designed to introduce a Beta- 
carbon at position 121 to displace the 6 -OH group of 

10 the third sugar residue from its hydrogen bonding 

position. In addition, the double mutant G121A/S411G 
was designed to investigate additivity of the two 
substrate specificity mutations. S411G is shown herein 
to reduce the ratio of initial rates of isomaltose 

15 production (from glucose condensation reactions) to 

that of glucose production (from the hydrolysis of DE 
10 maltodextrin) . 

The following provide further examples of the 
strategies used for the design of mutations having 

20 increased selectivity. 

300L.OQP mutation : According to the amino acid sequence 
homology study [Countinho and Reilly, 1994] , it was 
found that GAs from Rhizopus and some other fungal 
families have a longer amino acid sequence and form a 

25 larger loop or cavity in the S4 conserved region 

compared to A. niger or A. awamori GAs. Since single 
mutation events alone are unlikely to bring about 
substantial increase in the specificity of bond 
hydrolysis or synthesis, an insertional mutant GA was 

30 designed, designated 300Loop or 311-314Loop (SEQ ID 

NO: 2), and the inserted seven amino acids were adapted 
from Rhizopus GA because Rhizopus GA was the first 
enzyme to which the subsite theory was successfully 
applied [Himori et al., 1973]. The 300Loop mutation 

35 was designed to decrease the a- (1, 6) -glucosidic 

affinity by introducing a larger loop into the S4 



WO 98/03639 



PCTYUS97/12983 



-17- 

conserved region. 

Tyrl75Phe : Tyrl75 is within the third conserved region. 
The nearest distance between Tyrl75 and the fourth 
residue of inhibitor D-gluco-dihydroacarbose is 4.06 A 
5 [Stoffer et al. , 1995] . Tyrl75 is replaced by Phe or 
Gin in several other glucoaraylases . Changing Tyrl75 to 
Phe was designed to increase the hydrophobic 
interaction between enzyme and substrate. 
Glvl21Ala : Glyl21 is highly conserved in all 

10 glucoamylase sequences except in Clostridium sp. G005 
GA, which has high a- 1,6 activity and in which Gly is 
replaced by Thr. Since the <f> and \f/ of Gly 121 would 
allow an alanine in this position without causing a 
conformation distortion, G121A was designed to 

15 introduce a 0-carbon at position 121 to displace the 6- 
OH group of the third sugar residue from its hydrogen 
bonding position. 

Glyl21Ala with S411G (generally indicated as 
G121A/S411G) : The double mutant was designed to 

20 investigate additivity (cumulative) of the two 

substrate specificity mutations. S411G reduces the 
ratio of initial rates of isomaltose production (from 
glucose condensation reactions, see Examples) to that 
of glucose production (from the hydrolysis of 

25 maltodextrin 10) . 

The present invention provides a method of 
engineering mutations for fungal glucoamylase and then 
preparing engineered enzymes carrying cumulatively 
additive mutations. The initial step is to generate 

30 individual mutations by site directed mutagenesis and 
screen the individual mutations as described in the 
Examples. Those individual mutations which show 
decreased irreversible thermal inactivation rates or 
reduced isomaltose formation or increased pH optimum 

35 are then selected for combinational analysis. In 
general mutations are selected which have at least 
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wildtype reaction rates. 

Mutations are combined by site -directed 
mutagenesis to determine if their effects are additive 
as is discussed herein in the Examples. Site directed 
5 mutagenesis to produce enzymes carrying at least two of 
the isolated selected mutations is performed as is 
known in the art. These engineered enzymes are then 
screened for cumulatively additive effects on thermal 
stabilizing, pH optimum or reduced isomaltose 

10 formation. Alternatively the engineered enzymes 
carrying cumulative mutations are screened for 
cumulative effects on two or more of the parameters. 

For biochemical characterization of the mutants, 
GA is purified from culture supernatants of 15 -L batch 

15 fermentations by ultrafiltration, DEAE-Sephadex column 
chromatography, and column affinity chromatography 
using the potent inhibitor acarbose attached to a 
support [Sierks et al . , 1989]. Purities of the 
resulting preparations are tested by standard 

20 techniques such as SDS-polyacrylamide gel 

electrophoresis and isoelectric focusing with narrow- 
band ampholytes . Protein are measured by absorbance at 
280 nm or by Bradford's method [1976] . GA activity is 
measured by a glucose oxidase/o-dianisidine assay 

25 (Sigma) . . 

Selectivity is determined by any method known in 
the art but preferably by measuring the initial rate of 
isomaltose formation from 30% (w/v) glucose 
condensation reactions at pH 4.4 and 55°C in 0.05M 

30 sodium acetate buffer and then by measuring the initial 
rage of glucose formation in 30% (w/v) DE 10 
maltodextran hydrolysis reactions at pH 4.4 and 55°C 
0.05M sodium acetate buffer. From the resulting 
initial rates, the ratio of isomaltose formation to 

35 glucose formation is calculated. 

Thermostability is measured as is known in the art 
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but preferably by incubating the enzyme at selected 
temperatures between 65°C and 77.5°C at 2.5°C intervals 
followed by activity analysis at 35°C using 4% maltose 
as substrate. When first-order decay is observed, as 
5 with WT GA, decay rate coefficients are determined. 
Activation energies for decay are calculated from the 
rate coefficients at different temperatures. 

pH optimum is measured as is known in the art but 
preferably at 45°C at 16 pH values, ranging for 2.2 to 

10 7.0 using 0.025 M citrate -phosphate buffer with maltose 
or maltoheptaose as substrate. 

Saccharif ication is measured as described in the 
Examples. Briefly, glucoamylase is incubated with DE 
10 maltodextran as substrate in 0 . 05M sodium acetate 

15 buffer at pH 4.4 at 55°C. Samples are taken at various 
times from 0.5 to 288 hours and the production of 
glucose determined. 

The present invention provides vectors 
comprising an expression control sequence operatively 

20 linked to the nucleic acid sequence of the various 
mutant sequences disclosed herein, combinations of 
mutations and portions thereof. The present invention 
further provides host cells, selected from suitable 
eucaryotic and procaryotic cells, which are transformed 

25 with these vectors. 

Vectors can be constructed containing the cDNA of 
the present invention by those skilled in the art and 
should contain all expression elements necessary to 
achieve the desired transcription of the sequences. 

30 Other beneficial characteristics can also be contained 
within the vectors such as mechanisms for recovery of 
the nucleic acids in a different form. Examples are 
provided herein. Phagemids are a specific example of 
such beneficial vectors because they can be used either 

35 as plasmids or as bacteriophage vectors. Examples of 
other vectors include viruses such as bacteriophages, 
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baculoviruses and retroviruses, DNA viruses, cosmids, 
plasmids, liposomes and other recombination vectors. 
The vectors can also contain elements for use in either 
procaryotic or eucaryotic host systems. One of 
5 ordinary skill in the art will know which host systems 
are compatible with a particular vector. 

The vectors can be introduced into cells or 
tissues by any one of a variety of known methods within 
the art (calcium phosphate transf ection; 
10 electroporation; lipof ection; protoplast fusion; 

polybrene transf ection; ballistic DNA delivery; lithium 
acetate or CaCl transformation) . The host cell can be 
any eucaryotic and procaryotic cells, which can be 
transformed with the vector and which will support the 
15 production of the enzyme. 

The above discussion provides a factual basis for 
thermostable and selective mutants of fungal 
glucoamylase as well as methods of designing the 
mutations and screening for the cumulative effect of 
20 the mutations and vectors containing the mutations. 
The methods used with and the utility of the present 
invention can be shown by the following non- limiting 
examples and accompanying figures. 

EXAMPLES 

General methods in molecular biology: Standard 
molecular biology techniques known in the art and not 
specifically described were generally followed as in 
Sambrook et al. f Molecular Cloning: A Laboratory 
Manual, Cold Springs Harbor Laboratory, New York 
(1989), and in Ausubel et al., Current Protocols in 
Molecular Biology, John Wiley and Sons, Baltimore, 
Maryland (1989) and Rose, et al. Methods in Yeast 
Genetics: A Laboratory Course Manual, Cold Spring 
Harbor Laboratory Press, Cold Spring Harbor, NY (1990) . 
Polymerase chain reaction (PCR) was carried out 



25 



30 



35 
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generally as in PCR Protocols: A Guide To Methods And 
Applications, Academic Press, San Diego, CA (1990) . 
Oligonucleotides are synthesized as is known in the 
art. For example, an Applied Biosystems 380B DNA 
5 synthesizer can be used. 

Materials: 5. cerevisiae C468 (a leu2-3 leu 2-112 
his 3-11 his 3-15 mal") and the plasmid YEpPM18 were 
gifts from Cetus. Acarbose was a gift from Miles 
Laboratories. All restriction enzymes were purchased 

10 from Promega as well as T4 DNA ligase and pGEM-7Z(+), 

an E. coli phagemid vector, were from Promega. Maltose 
(G 2 ) , maltotriose (G 3 ) , maltotetraose (G,) , raalto- 
pentaose (G 5 ) , maltohexaose (G 6 ) , maltoheptaose (G 7 ) , 
glucose oxidase, peroxidase, and a-naphthol were from 

15 Sigma. Isomaltose (iG 2 ) was purchased from TCI America. 
DE 10 Maltodextrin with the average degree of 
polymerizations (DP) of 10, 6, and 4, respectively, 
were from Grain Processing Corporation. High- 
performance thin- layer chromatographic (HPTLC) plates 

20 (LHPK silica gel €0 A, 20 x 10 cm) were from Whatman. 
Site-directed mutagenesis: Site-directed 
mutagenesis was performed according to the Muta-Gene 
phagemid in vitro mutagenesis kit from Bio-Rad which is 
based on the method of Kunkel et al [1985] . A 1.7 kb 

25 XhoI-»BamHI DNA fragment coding for the glucoamylase 

catalytic domain was cloned into a pBluescript II KS(+) 
vector from Stratagene. Oligonucleotides used as 
mutagenic primers are provided with the specific 
Example. The presence of the individual mutations was 

30 confirmed by sequencing and each mutated GA gene 

fragment was subcloned into YepPM18 [Cole, et al., 
1988] and transformed into S. cerevisiae. 

Enzyme production and purification: Wild- type (WT) 
and mutant enzymes are produced by growing yeast at 

35 30°C in 5.3 L SD + His media for 72 hours at pH 4.5 in 
a 5.0 L fermentor. After 48 hours, lOOg of dextrose 



WO 98/03639 



PCMJS97/12983 



-22- 

and 22g of (NH 4 ) 2 S0 4 in 300ml H 2 0 is added as a 
supplement. Following growth, the culture is 
centrifuged to remove yeast cells, the supernatant is 
concentrated by ultrafiltration, diafiltered against 
5 0.5 M NaCl/0.1 M NaOAc , pH 4.5 and purified by 

acarbose-sepharose affinity chromatography. GA is 
eluted with 1.7 M Tris-Cl, pH 7.6, dialyzed against H 2 0, 
further concentrated by ultrafiltration and diafiltered 
against 0.05 M NaOAc buffer, pH 4.5. The protein 

10 concentration is determined according to the Pierce 
bicinchoninic acid protein assay [Smith et al., 1985] 
using bovine serum albumin as a standard. 

Enzyme activity assays: Enzyme activities were 
determined at 50°C using 4% maltose in 0.05 M NaOAc 

15 buffer pH 4.5 as substrate. One international unit 
(IU) of enzyme activity was defined as the amount of 
enzyme required to produce 1 /imol/min glucose at assay 
conditions. Following mixing enzyme with substrate, 
six 100 /il samples were removed at seven minute 

20 intervals over 42 minutes, the reaction stopped with 40 
pi of 4.0 M Tris-Cl, pH 7.0 and the glucose 
concentration was determined according to the Sigma 
peroxidase -glucose oxidase/o dianisidine glucose assay 
kit. 

25 Irreversible thermal inactivation : Duplicate 

aliquotes of 40 fig/ml of purified wild- type and mutant 
enzymes were subjected to inactivation at six or more 
temperatures between 65° and 80°C at intervals of 
2.5°C. Samples were removed at six different time 

30 points, immediately placed on ice and stored at 4°C for 
24 hours. The residual activity of the inactivated 
samples along with a corresponding sample which had not 
been subjected to thermal inactivation, was determined 
as described above but at 35°C. 

35 pH dependence of glucoamylase activity: pH 

dependence of glucoamylase activity was measured at 
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45°C at 16 different pH values, ranging from 2.2 to 
7.0, using 0.025 M citrate -phosphate buffer [Mcllvane, 
1921] with maltose or maltoheptaose as substrate. The 
ionic strength of the citrate-phosphate buffer was 
5 maintained at 0.1 by adding potassium chloride. The pK 
values of free enzyme and enzyme -substrate complex were 
measured at substrate concentrations (i) smaller than 
0.2 ^, so that the initial rate (v) was proportional to 
^cat/X/ and (ii) higher than 10 K„, so that the initial 

10 rate (v) was proportional to k cat [Sierks & Svensson, 

1994, see also Whitaker (1994) Principle of enzymology 
for the food sciences, 2 nd edition. Marcel Dekker, NY] . 
The pK values of two catalytic groups of free enzyme 
and enzyme- substrate complex were calculated by fitting 

15 the initial rates as a function of pH values to the 
equation log Y - log [C/ (1+H/Kj+K 2 /H] by using the 
software of Enzfitter. Y is the observed value of the 
parameter of interest (i.e. k cat /K m or k cat ) measured at 
different pH values, C is the pH independent value of Y 

20 (i.e. the maximal value of k cat /K m or k cac ) , H is the 
concentration of hydrogen ion, K x and K 2 are 
dissociation constants of catalytic groups of enzyme. 
When the values of apparent pK x and pK 2 were separated 
by less than 3 pH units, the pK values were adjusted by 

25 equations + (JT) 2 - + 4(fT) opt and (H*) opt = <l K X K 2 

[Whitaker, 1994] . The concentration of hydrogen ion at 
the optimum pH, (H*)^, was calculated from pHopt which 
is equal to the average of apparent pK x and pK 2 . (JT) . x 
and (fT) 2 (apparent and K 2 ) correspond to the 

30 concentrations of hydrogen ion when the pH values are 
equal to apparent pK x and pK 2 , respectively. 

The hydrolysis of DK 10 maltodextrin 
(Saccharification) : Hydrolysis was performed at 35°C 
and/or 55°C (as indicated in the text) with 28% (w/v) 

35 DE 10 maltodextrin as substrate in 0.05 M sodium 

acetate buffer at pH 4.4 with the addition of 0.02% 
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sodium azide, used to inhibit microbial growth in the 
reaction mixtures. The enzyme concentration was 2.64 
/iM for both wild-type and mutant GAs . Samples were 
taken at various times (from 0.5 to 288 hours) and the 
5 reactions were stopped by adding samples to the same 
volume of 1 M Tris-HCl buffer at pH 7.0, since Tris is 
a known inhibitor of glucoamylase [Clarke & Svensson, 
1984] . The production of glucose was determined by the 
glucose oxidase method [Rabbo & Terkildsen, 1960] . 

10 Initial rates of glucose production were determined by 
fitting the experimental data to the equation c = 
At/ (1+Bt) , where c is the product concentration, t is 
time, and A (the initial rate) and B are obtained from 
the nonlinear regression. At 55°C, only the time 

15 points before 70 hours were used for the calculations, 
since the glucose production by that time had already 
declined for wild-type GA. 

Glucose condensations reactions: Glucose 
condensation reactions were performed at 35°C and 55°C 

20 with 30% (w/v) D-glucose as substrate in 0.05 M acetate 
buffer at pH 4.4 for 12 days with the addition of 0.02% 
sodium azide, used to inhibit microbial growth in the 
reaction mixtures. The enzyme concentration was 2.64 
fM for both wild-type and mutant GAs . Samples were 

25 taken at various times and the reactions were stopped 
by adding samples to the same volume of 1 M Tris-HCl 
buffer at pH 7.0. High Performance Thin Layer 
Chromatography (HPTLC) and Imaging Densitometry were 
used to determine the production of isornaltose by a 

3 0 method modified from that described by Robyt et al. 
[Robyt and Mukerjea, 1994] . One microliter of 
variously diluted samples and six different 
concentrations of standard (containing glucose, maltose 
and isornaltose) were applied to the HPTLC plates. The 

35 developing solvent system contained acetonitrile, ethyl 
acetate, 1-propanol and water in the volume proportions 
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of 85:20:50:40. Only one ascent was used to develop 
the carbohydrate separation on HPTLC plates. After 
development, the plates were air-dried, dipped into an 
EtOH solution containing 0.3% (w/v) alpha-naphthol and 
5 5% (v/v) H 2 SO«, air-dried again, and incubated 
approximately 10 min at 120°C to visualize the 
carbohydrates. Densities of the isomaltose spots on 
HPTLC plates were quantified by Imaging Densitometry 
(Bio-Rad, Model GS-670) , using Molecular Analyst 
10 software (Bio-Rad) . The experimental data were fitted 
to the equation c = At/(1+Bt), described above for the 
hydrolysis of DE 10 maltodextrin, to obtain the initial 
rates of isomaltose production. 

15 

EXAMPLE 1 

STABILIZATION OF ASPERGILLUS AWAMORI GLUCOAMYLASE 
BY PROLINE SUBSTITUTION 

20 The following example is an exemplar of the 

methods and procedures that are used in the analysis of 
an individual mutation of a glucoamylase . To 
investigate the mechanisms governing Aspergillus 
awamori glucoamylase thermal stability, three proline 

25 substitution mutations were constructed. These 

mutations were predicted to increase GA stability by 
decreasing the enzyme's conformational entropy of 
unfolding . 

Aspergillus awamori glucoamylase (or-1, 4-D-glucan 
30 glucohydrolase, EC 3.2.1.3; GA) is an enzyme which 
catalyses the release of /3-glucose from the non- 
reducing ends of starch and related oligosaccharides. 
GA is used in, and defines the rate limiting step of, 
the commercial conversion of starch to high glucose 
35 syrups which may be converted to fructose syrups by 

glucose isomerase, or used in fermentations to produce 
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ethanol. GA is used industrially at 55°-60°C; at 
higher temperatures the enzyme is rapidly and 
irreversible inactivated. Therefore, a GA 
variant with increased thermostability would be 
5 advantageous industrially to decrease reaction times 
and/or to increase solids concentrations. 

Previous work has shown that the natural stability 
of oligo 1,6-glucosidase [Suzuki et al., 1987] and 
pullulanase [Suzuki et al., 1991] can be positively 

10 correlated to the mole percent proline present in the 
protein, and a general rule for protein stability has 
been proposed [Suzuki, 1989] . This work has been 
extended to show that bacteriophage T4 lysozyme 
[Matthews et al., 1987] and Bacillus cereus ATCC 7064 

15 oligo 1,6 glucosidase [Watanabe et al, 1994] can be 

stabilized by engineering proline into selected sites 
thereby decreasing the protein's conformational entropy 
of unfolding. 

Three sites (Ser30, Asp345 and Glu408->Pro) were 

20 selected for proline substitution based on structural 
and evolutionary considerations. Mutations at these 
sites were constructed using the cloned A. awamori gene 
[Innis et al, 1985] and the proteins were expressed in 
Saccharomyces cerevisiae [Cole, et al, 1988] . The 

25 stability of the mutant proteins was measured by their 
resistance to irreversible thermal inactivation at 
various temperatures. As shown herein, the Ser30-*Pro 
mutation increased. However, unexpectedly the 
Glu408-*Pro mutation decreased and the Asp345-*Pro 

30 mutation did not significantly alter GA stability. 
Site-directed mutagenesis: Site-directed 
mutagenesis was performed as described herein above. 
The following oligonucleotides were used as mutagenic 
primers: CAGAGTCCGCGCCCGGCACCCAAGCACCGTC (Ser30-*Pro) 

35 (SEQ ID No: 3) , AAGTCCAGCGACACAGGTGTGACCTCCAACGAC 
(Asp345-»Pro) (SEQ ID No: 4) and CGAGCGGAAAGCTGC 
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GGGCCATCAGACTTGTC (Glu408->Pro) (SEQ ID No: 5) . 

Selection of sites for proline substitution: 

Based on the nearly identical catalytic domain of A. 
awamori var X100 GA whose structure is known [Aleshin 
5 et al # 1992] three sites for substitution were chosen, 
which met the following criteria: 1) Ramachandran 
(0,^) angles were within allowed values for proline 
[Ramachandran et al., 1963] . For this work the <f> and \p 
angles at the substituted site were restrained to the 

10 broad range 4> = -90° to -40°, ^ = 120° to 180° or 

<f> = -90° to -40° , \p » -50° to 10°. 2) Residues were 
highly solvent exposed, since mutation of residues in 
the core of the enzyme were thought to be more likely 
to decrease the enzyme's catalytic efficiency. 3) 

15 Residues didn't participate in hydrogen bonding with 
other amino acids. Additionally, based on sequence 
alignments with GA's from other organisms [Coutinho and 
Reilley, 1994b] only residues which met the above 
structural criteria and were hot well conserved were 

20 selected for mutation. Ser30 could be aligned with 
proline in GAs from Hormoconis grisea var thermoidea 
and H. resiae GamP [Coutinho and Reilly, 1994b] , which 
made it particularly attractive for proline 
substitution. 

25 

RESULTS 

Specific activity 

None of the proline substitution mutations 

significantly altered enzyme specific activities of 
30 wild-type and mutant GA's at 50°C and pH 4.5. This 

suggests that these mutations did not significantly 

alter the enzyme's structure around the active site or 

alter its interaction with substrate. 

Irreversible thermal stability 
35 Wild-type and mutant GA's were subjected to 

thermal inactivation at pH 4 . 5 as described in the 
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experimental protocols. Serai logarithmic plotting of 
the percent residual activity against inactivation time 
yielded inactivation rate coefficients {kd) . Figure 1 
shows the relationship between temperature and Jed for 
5 wild-type and mutant GA's. Based on these data, 

activation energies for thermal inactivation (AG 1 ) were 
calculated using transition state theory and melting 
temperatures (Tm) , the temperature at which the enzyme 
was 50% inactivated after 10 minutes were computed 

10 (Table 1) , As can be seen, the Glu408-*Pro mutation 
greatly decreased, the Asp345->Pro mutation did not 
significantly alter and the Ser30-*Pro mutation 
increased GA stability . 

It should be noted that although Table 1 shows 

15 that the Asp345-*Pro mutant GA demonstrated slightly 
increased AG 4 and Tm, these changes are generally not 
significant or that the Asp345-»Pro mutant GA is more 
stable than wild- type since the Jcds for this mutant 
enzyme at two well separated temperatures (65° and 

20 75 °C) are essentially indistinguishable from wild- type 
(Figure 1) . 

The proline substitution mutations had different 
thermostabilities when measured by their resistance to 
irreversible thermal inactivation. When compared to 

25 wild-type GA, Glu408-*Pro decreased, Asp345-*Pro did not 
significantly alter and Ser30->Pro increased GA 
stability (Figure 1 and Table 1) . 

Glu408-*Pro destabilized GA. As was first 
suggested by Schimmel and Flory [1968] and has been 

30 expanded by others [MacArthur and Thornton, 1991; 

Hurley et al, 1992] proline not only restricts the 0,^ 
values for the site at which it exists, but also the 
0,^ values of the preceding residue. These reports 
suggest that the (<f>,\l>) values for the residue preceding 

35 proline should be restricted to approximately 0 » -180° 
to -55° and 0= 55° to 180° or 0 = -180° to -55° and 
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\f/= -30° to -70° for all residues in Xaa-Pro except for 
Xaa r Gly, for which the preceding still applies, but is 
extended to include #=45° to 180°. In the published A. 
awamori var. X100 catalytic domain structure [Aleshin 
5 et al., 1992], Asp408 (#=-65°, ^=146°) which aligns 
with Glu408 in A. awamori GA, has <f>,)// values within 
ranges acceptable for proline. However, the preceding 
residue Gly407 (#=80°, ^=-5°) has outside 
acceptable ranges for positions preceding proline. It 

10 is not surprising then, that the Glu408-»Pro 

destabilized GA. Additionally, X-ray crystallography 
suggests that position 408, in the closely related A. 
awamori var. X100 GA 2 , lies within a /?- strand; a site 
not well suited for proline substitution. 

15 Asp345 (<f> = -65°, ^ = -26°) and the preceding 

Thr344 (0 - -116°, \p =178°) have <f>,\{/ angle values lay 
well within allowed values for proline substitution at 
position 345. However, the Asp345->Pro mutant GA did 
not demonstrate stability significantly different from 

20 wild- type GA. This is particularly unexpected since 
position 345 lies at the N-terminus of an a-helix 2 ; a 
position previously shown to be particularly favorable 
for proline substitution [Watanabe et al, 1994] . 

Ser30 (<f> = -49°, \p= 130°) is preceded by Val29 {0 

25 = -127° , \p = 46°) both of which have acceptable 4>,^ 
angle values except Val29 V = 46° which is slightly 
smaller than ideal for proline substitution at position 
30. 

In summary, when expressed in Saccharowyces 
30 cerevisiae, Glu408-*Pro greatly decreased, Asp345-»Pro, 
did not significantly alter and Ser30-»Pro strongly 
stabilized the enzyme. The Ser30-»Pro mutant GA showed 
a significantly decreased rate of irreversible thermal 
inactivation when analyzed between 65° and 77.5°C 
35 without decreased enzyme activity. At 65 °C a 1.7- fold 
decrease in thermal inactivation rate coefficients was 
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seen and the activation energy for thermal inactivation 
was increased by 1.6 kJ/mol relative to wild- type GA. 

EXAMPLE 2 

5 ENGINEERED DISULFIDE BONDS 

The following example is an exemplar of the 
methods and procedures that are used in the analysis of 
an individual mutation of a glucoamylase „ The process 
of GA thermoinactivation is thought to be dominated by 

10 formation of enzymes with incorrect conformation [Munch 
and Tritsch, 1990] . Previous work supported this 
hypothesis. Site-directed mutagenesis has been used to 
eliminate sites of deamidation and peptide hydrolysis 
Chen et al . , 1994 a,b) . The corresponding mutations 

15 Asnl82-*Ala and Asp257-»Glu had reduced irreversible 
thermoinactivation rates at pH 4.5 below 70°C but 
increased rates above 70°C. Thus GA thermoinactivation 
is predominantly caused by "scrambled" structures 
rather than by deamidation and peptide hydrolysis. 

20 Furthermore, mutations Glyl37-*Ala, Glyl3 9->Ala and 

Glyl37/l39-*Ala/Ala, made to reduce helix flexibility, 
showed increased thermostability up to 75°C (Chen et 
al., 1996) apparently by slowing down the formation. of 
incorrect structures. 

25 To improve protein thermostability by preventing 

formation of incorrect structures, several strategies 
have been proposed including introducing covalent 
linkage such as disulfide bonds (Perry and Wetzel, 
1984; Wetzel, 1987; Matsumura et al., 1989, Clarke and 

30 Fersht, 1993) . 

There are a total of nine cysteine residues in A. 
awamori GA, eight of which form disulfide -linked pairs, 
which are assumed to enhance the folding and stability 
of GA, residues 210 and 213, 262 and 270, 222 and 449 

35 [Aleshin et al., 1992] and 509 and 604 [Williamson et 
al. f 1992b]. In this Example, additional disulfide 
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bonds are introduced into GA to explore the effect on 
thermostability and catalytic activity. Two engineered 
disulfide bond mutants designated A27C/N20C 
(abbreviated S-S) and A471C/T72C were constructed. The 
5 new disulfide bond formed by A27C/N20C connects the C- 
terminus of helix 1 (Asn20) and a turn where residue 
Ala27 is located, while A471C/T72C bridges the N- 
terminus of helix 3 and the end of the 30 -residue 
highly O-glycosylated belt region together. The 

10 disulfide bonds are formed spontaneously after 
fermentation and have different effects on GA 
thermostability and catalytic activity. v 
Site-Directed Mutagenesis: Site-directed mutagenesis 
was performed as described herein above. 

15 Oligonucleotide primers used are: 5'-CGT ACT GCC ATC 
CTG TOT AAC ATC GGG GCG GA-3' (N20C, AAT-*TGT) (SEQ ID 
No: 6) , 5 ' -ATC GGG GCG GAC GGT TGT TGG GTG TCG GGC GCG- 
3' (A27C, GCT-VTGT) (SEQ ID No: 7), 5'-CGA AAT GGA GAT 
TGC AGT CTC-3' (T72C, ACQVTGC) (SEQ ID No:8), 5 f -G AGT 

20 ATC GTG TGT ACT GGC GGC ACC-3' (A471C, GCT-VTGT) (SEQ ID 
No: 9), with the underlined letters indicating the 
nucleotide mutations. 

SDS-PAGE and Thio- titration: SDS-PAGE was carried out 
using 0.75 mm thick 10% polyacrylamide gels following 

25 the method of Garfin [1990] . For thio- titration, GA at 
2 mg/ml concentration was denatured by boiling in 
denaturing solution containing 2% SDS, 0.08 M sodium 
phosphate (pH 8.0) and 0.5 mg/ml EDTA [Habeeb, 1972] 
with or without 50 mM DTT [Pollitt and Zalkin, 1983] 

30 for 10 min. The denatured GA (reduced or non-reduced) 
was concentrated using Centricon 30 concentrators 
(Ami con, MA, USA) and the reduced GA was applied to 
Bio-spin 30 chromatography columns (Bio-Rad, CA, USA) 
pre-equilibrated with denaturing solution to remove 

35 DTT. The resulting solution as well as the non-reduced 
denatured GA sample were divided into two portions. 
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One portion was used for a protein concentration assay 
and the other portion was assayed for thio reduction by 
mixing with 4 mg/ml DTNB in denaturing solution with a 
30:1 volume ratio, followed by incubation at room 
5 temperature for 15 minutes, and absorbance measurement 
at 412 nm with a molar absorptive value of 13,600 M^cm" 1 
[Habeeb, 1972] . 

GA Activity Assay: As described herein above, 
maltose was used as substrate in enzyme kinetics 

10 studies, with concentrations ranging from 0.2 to 4 Kp 
at 35°C and pH 4.5 as described previously [Chen et al . , 
1994b] . Kinetics parameters were analyzed by the 
program ENZFITTER. In residual enzyme activity assays, 
the conditions are the same as in the enzyme kinetics 

15 studies except that only one concentration of maltose 
(4%) is used as substrate. Specific activity assays 
were carried out with 4% maltose as substrate at 50°C 
and pH 4.5. One unit (IU) was defined as the amount of 
enzyme required to produce 1/imol glucose per min under 

20 the conditions of the assay. To compare the 

temperature optima of catalytic activities of wild-type 
and mutant GA, activities were assayed at pH 4.5 with 
4% maltose as substrate at different temperatures. 
Irreversible Thermo inacti vat ion : As described 

25 herein above, purified wild- type or mutant GA proteins 
were incubated at five different temperatures from 65°C 
to 75°C at 2.5°C intervals at 40 /xg/ml in 0.05 M NaOAC 
buffer (pH 4.5). At six different time points, 
aliquots of the incubating enzyme were removed, quickly 

30 chilled on ice, stored at 4°C for 24 hours, and 

subjected to residual activity assay. The irreversible 
thermoinactivation of GA obeyed first-order kinetics 
[Chen et al., 1994b]. Thermoinactivation rate 
coefficients, k d were determined as described previously 

35 [Chen et al., 1994b]. 

Computer Modeling and Three-dimensional View of 
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Mutated Residues: The candidate residues of A. 
awamori GA to form disulfide bond were modeled with the 
crystal structure of A. awainori var. X100 GA [Aleshin 
et al., 1992] (Igly in the Brookhaven Protein Data 
5 Bank) as reference by the SSBOND program (Hazes and 
Dijkastra, 1988) installed in a DEC 3100 workstation. 

Selection of Mutation Site: Residues Asn20, Ala27 
and Thr72, Ala471 were chosen to be replaced with 
cysteine. After the analysis of crystal structure of 

10 A. awamori var. X100 GA [Aleshin et a 1 . , 1992] by the 
program SSBOND, 132 pairs of residues were found that 
could potentially be sites for a disulfide bond. Pairs 
containing glycine were discarded on the assumption 
that glycine may be required for flexibility at that 

15 site. Also, the residues involved in hydrogen bonds 
and electrostatic interactions were eliminated. 
Residues 20 paired with 27 as well as 72 paired with 
471 were chosen as candidates for disulfide bond 
formation according to the geometrical analysis . Amino 

20 acid sequence alignment among related GAs showed that 

there is a disulfide bond between position 20 and 27 in 
Neurospora crassa [Coutinho and Reilly, 1994b] , which 
suggested that introducing disulfide bond between 
position 20 and 27 would not cause unfavored 

25 interactions there in A. awamori GA. Furthermore, the 
20/27 disulfide bond would link the C-terminus of helix 
1 and the conserved SI fragment of GA involved in 
substrate binding [Coutinho and Reilly, 1994a] to form 
a loop, near another loop very critical for catalysis 

30 containing Trp 120, a residue involved in substrate 

binding [Sierks et al., 1989]. Therefore, the proposed 
20/27 disulfide bond was expected to stabilize GA by 
keeping the correct conformation for catalysis and 
substrate binding. 

35 Another further candidate for a disulfide bond 

pair was between positions 471 and 72. This disulfide 
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bond would link the N- terminus of helix 3 and the end 
of the 30 -residue (440-470) highly O-glycosylated belt 
region to form a loop. This disulfide bond also would 
make an additional linkage between the catalytic domain 
5 and the O-glycosylated linker. This O-glycosylated 
linker has been proved to be important for GA 
thermostability by limiting the conformational space 
available to the GA unfolded peptide [Semimaru et al. # 
1995 and Williamson et al., 1992]. This disulfide bond 

10 could have a globe effect on the thermostability of GA 
because of this linkage. The side chain -OH group of 
Thr72 in A. awamori var. X100 GA [Aleshin et al., 1992] 
is hydrogen bonded to the main chain N atom of Asp73 . 
In A. awamori GA, however serine is found at residue 73 

15 in place of Asp. It is possible that the hydrogen bond 
between residues 72 and 73 does not exist in A. awamori 
GA, and therefore replacing Thr72 with Cys will not 
disturb this interaction. This hydrogen bond is 
apparently not critical for GA since Thr72 is replaced 

20 by Ala, Lys or Val in other GAs [Coutinho and Reilly, 
1994b] . 

Engineered Disulfide Bonds Were Formed 
Spontaneously: After GA purification, the engineered 
disulfide bonds were found to be formed spontaneously 

25 by the following two approaches. 

First, the mutant A471C/T72C has faster mobility 
than wild-type during SDS-PAGE under non-reducing 
, conditions, suggesting that an additional disulfide 
bond forms a new loop retarding the migration. The 

30 possibility that a truncated enzyme was formed in this 
case was eliminated by DNA sequencing of the mutant 
cDNA and MALDI analysis. The MALDI data showed that 
the mutant GA had the same molecular weight as wild- 
type GA. Mutant A27C/N20C has the same migration as 

35 wild- type GA, which may be because the additional loop 
caused by the engineered disulfide bond is too small 
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( seven residues) to affect migration. 

Second, the new disulfide bonds were demonstrated 
by thio group titration. Comparing the numbers of free 
thio groups before and after the treatment of reducing 
5 reagent DTT, the total disulfide bonds in mutant and 
wild-type GA were deduced as reported in Table 2. 
Wild- type, A27C/N20C, and A417C/T72C GA have in total 
8.6, 10.9, and 10.4 free thio groups respectively 
according to the [SH] /molecule ratio in the presence of 

10 reducing reagent DTT (Table 2). In the absence of DTT 
the numbers are 0.9, 0.9 and 1.3, respectively (Table 
2) . This suggested that the number of disulfide bonds 
among wild- type, A27C/N20C and A471C/T72C are 4, 5 and 
5, respectively. Therefore, the introduced cysteine 

15 residues formed disulfide bonds instead of remaining 
free thiols. 

Enzymatic Activity and Optima Temperature of 
Catalysis: The enzymatic properties of double 
mutations A27C/N20C and A471C/T72C were not changed 

20 compared with wild-type at 35°C and 50°C as shown in 

Table 3, while single mutations had significant reduced 
activity. Mutant A27C/N20C and A471C/T72C had specific 
activities at 50°C and kinetic parameters at 3 5°C very 
close to wild- type GA (Table 3) . The single mutant 

25 A27C had slightly increased K,, but the same Jc catc value 
as wild-type GA, and thus a reduced k cat /K m ratio of 
-30%. Mutant N20C had the same K m but both a decreased 
k cat and k cat /K n ratio and a decreased specific activity 
at 50°C of more than 50%, 

30 Irreversible Thermoinactivation of GA: The 

irreversible thermoinactivation of wild- type and mutant 
GA was studied at 65°C, 67.5°C r 70°C, 72.5°C and 77.5°C 
with first-order irreversible theremoinactivation 
coefficients k d shown in Figure 2. Mutants A27C, 

35 A27C/N20C and A471C/T72C had smaller k d values than did 
wild- type GA within the measured temperature range, 
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which means the activity decayed more slowly than wild 
type, while mutant N20C had greater k d value than wild- 
type at all temperatures except 75°C, which means that 
N20C decayed faster than wild-type. 
5 Table 4 shows the activation enthalpy (AH t ) , 

entropy (AS t ) and free energy of unfolding (AG,) at 65°C 
and 75°C of wild-type and mutant GAs, calculated 
according to transition- state theory. The enthalpies 
of N20C and A27C/N20C decreased by 42 and 24 KJ/mol 

10 respectively, while no significant change occurs for 
A27C and A471C/T72C. Mutants N20C and A27C/N20C had 
decreased entropy of 115 kJ/mol and 75 kJ/mol 
respectively, while entropy of mutants A27C and 
A471C/T72C showed no significant change. Mutant A27C 

15 and A471C/T72C had a slightly higher AG* than wild- type 
GA at 65°C and 75°C (<0.5 kJ/mol) , while the AG 1 of 
A27C/N20C was higher than that of wild-type by 1.5 and 
2.2 kJ/mol at 65°C and 75°C respectively. Mutant N20C 
had a decreased AG* by 3 . 0 and 1.8 kJ/mol at 65°C and 

20 75°C, respectively, compared with wild-type GA. 
Therefore, the engineered disulfide bond mutant 
A27C/N20C significantly increased GA thermostability 
compared with wild- type GA while the single mutants 
produced either a slight increase (A27C) or a slight 

25 decrease (N20C) in thermostability. The other 

disulfide bond mutant had the thermostability identical 
to wild- type GA. 

EXAMPLE 3 

MUTATION A27C/N20C IN COMBINATION WITH OTHER MUTATIONS 

30 

In previous studies applicants have constructed 
the thermostable mutants G137A [Chen et al., 1996] and 
S436P (Li et al., 1996), which have the potential to be 
combined and improve thermostability additively. In 
35 this Example, these mutations are combined with each 

other and with A27C/N20C (S-S; Example 2) to test their 
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effects (cumulative/additive) on thermostability and GA 
activity. 

Enzymatic Activity and Optima Temperature of 
Catalysis: The combined mutants A27C/N20C/G137A and 
5 A27C/N20C/S436P had increased specific activity while 
mutant G137A/S436P had similar specific activity to 
wild-type GA (Table 3) . The double mutants A27C/N20C 
and A471C/T72C as well as the combined mutant 
A27C/N20C/G137 had changed optimal temperatures for 

10 catalysis. 

Relative activity assays at temperatures from 60°C 
to 74°C (Figure 3) showed that wild-type, mutant 
A27C/N20C and A471C/T72C had the highest activity at 
71°C, 72°C and 72.5°C, respectively. From 60°C to 

15 67.5°C, mutant and wild-type GA had very similar 

activities. However, when the temperature was above 
70°C, their relative activities differed substantially. 
Mutants A27C/N20C and A27C/N20C/G137A had higher 
activity than wild-type consistently from 70°C to 76°C 

20 with a peak at 72.5°C, while mutant A471C/T72C had 

activity lower than wild-type from 70°C to 71°C and 73°C 
to 74°C but higher at 72°C which is its optimal 
temperature. Thus mutant GAs A27C/N20C, A471C/T72C and 
the combined mutant A27C/N20C/G137A had increased 

25 temperature optima above wild- type GA by 1.5°C. 

Irreversible Thermoinactivation of GA: The 
irreversible thermoinactivation of wild- type and mutant 
GA was studied at 65°C, 67.5°C, 70°C / 72.5°C and 77.5°C 
with first-order irreversible theremoinactivation 

30 coefficients k d shown in Figure 2. Mutants A27C, 
A27C/N20C and A471C/T72C, A27C/N20C/G137A, 
A27C/N20C/S436P and G137A/S436P had smaller k d values 
than did wild-type GA within the measured temperature 
range, which means the activity decayed more slowly 

35 than wild type, while mutant N20C had greater k d value 
than wild-type at all temperatures except 75°C, which 
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means that N20C decayed faster than wild- type. 

Table 4 shows the activation enthalpy (AH>) , 
entropy (AS t ) and free energy of unfolding (AG,) at 65°C 
and 75°C of wild- type and mutant GAs, calculated 
5 according to transition-state theory. 

The helix flexibility mutant Gl 3 7A showed additive 
thermostability when combined with either S436P or 
A27C/N20C. The combination S436P with A27C/N20C did 
not show additivity. 

10 

EXAMPLE 4 

FURTHER STUDIES WITH COMBINED MUTATIONS 

To further investigate whether individual 
stabilizing mutations can cumulatively stabilize 

15 Aspergillus awamori glucoamylase (GA) , mutant enzymes 
were constructed containing combinations of 
thermostabilizing mutations. Previous work has shown 
that the following mutations stabilize GA as 
demonstrated by decreased irreversible thermal 

20 inactivation rates when inactivated in the absence of 
carbohydrate: Ser30-»Pro (S30P; Example 1), Glyl37->Ala 
(G137A), and Asn20-Cys/Ala27->Cys (which creates a 
disulfide bond between residues 20 and 27 and is 
therefore noted as S-S for convenience; Example 2) . To 

25 investigate whether individual stabilizing mutations 
can cumulatively stabilize GA # additional combined 
mutant enzymes were prepared utilizing the these three 
mutations. 

Site-directed mutagenesis: The S-S/S30P/G137A combined 
30 mutant was constructed using the S-S/S30P 

oligonucleotide listed above and a single stranded DNA 
template derived from a pBluescript II KS(+) vector 
with a 1.7 Jcb XhoI-*BamHI DNA fragment coding for the GA 
catalytic domain which already contained mutations 
35 conferring the S30P and G137A amino acid substitutions. 
The presence of the individual mutations was confirmed 
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by sequencing and each mutated GA gene fragment was 
subcloned into YEpPM18 [Cole et al., 1988] and 
transformed into S. cerevisiae. 

Thiol analysis: 10 nmol of wild- type, S-S/S30P and S- 
5 S/S30P/G137A mutant GAs were incubated in 0.2 mM 5,5'- 
dithiobis (2-nitrobenzoic acid), 6 M GdnHCl, and 50mM 
Tris, pH 8 in duplicate [Fierobe et al., 1996]. The 
thiol concentration was calculated from a standard 
curve established using 0-30/xM cysteine. 

10 Irreversible thermal inactivation : Wild- type and 

mutant GAs were subjected to thermal inactivation at 
six or seven temperatures between 65° and 80°C at 2.5°C 
intervals in duplicate. Following 24 hours at 4°C, the 
residual activities of the inactivated samples were 

15 analyzed at 35°C along with a corresponding sample 
which had not been inactivated [Chen et al, 1996] . 
Saccharif ication analysis: Saccharif ications were 
performed in duplicate using stirring heating blocks 
and tightly sealed vials to prevent evaporation. Eight 

20 fig/ml of wild-type and mutant GAs were assayed using 

28% (w/v) Maltrin DE 10 maltodextrin in 0.05 M NaOAc pH 
4.5 as substrate. At various times, sample was 
removed, diluted appropriately in 0.05 M NaOAc pH 4.5 
and the reaction stopped by adding 100/xl of diluted 

25 sample to 40 /il of 4.0 M Tris-Cl, pH 7.0. The glucose 
concentration was determined by a glucose oxidase/o 
dianisidine assay [Banks and Greenwood, 1971] . 

RESULTS 
Enzyme activities 

30 Table 5 shows the specific activities of the wild- 

type and mutant GAs at 50 °C and pH 4.5 using maltose as 
substrate. None of the mutant GAs demonstrated reduced 
enzyme activity and the S30P/G137A and S-S/S30P/G137A 
mutants were somewhat more active than wild-type at 

35 50°C. To further investigate this observation, the 
activities of these mutant enzymes were assayed at 
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various temperatures between 35° and 68°C {Figure 4) . 
The S30P/G137A and S-S/S30P/G137A mutant OAs were more 
active than wild-type at all temperatures examined. 
Thiol analysis 
5 The formation of a disulfide bond between 

positions 20 and 27 in the Asn20->Cys/Ala27-*Cys mutant 
GA has been confirmed (Example 2) . Table 6 shows the 
results of thiol analysis for the combined mutants S- 
S/S30P and S-S/S30P/G137A. A. awamori GA has one free 

10 cysteine at position 320. The combined mutant GAs show 
slightly higher thiol content per molecule than wild- 
type which may reflect less than complete disulfide 
bond formation between positions 20 and 27. However, 
if the disulfide bond was completely unformed, the 

15 [SH] /protein would be expected to rise to approximately 
three with the addition of two free cysteine residues. 
Therefore, we conclude that the disulfide bridge is 
formed at 70-80% the expected theoretical yield for 
complete formation. 

20 Irreversible thermal inactivation 

Wild- type and mutant GAs were subjected to thermal 
inactivation at pH 4.5 between 65° and 80°C. 
Semilogarithmic plotting of residual activity versus 
inactivation time yielded inactivation rate 

25 coefficients {kd) . Figure 5 shows the effect of 

temperature on kd for wild-type and mutant GAs. As can 
be seen, the combined mutants are significantly more 
stable than the individual mutant enzymes. 
Additionally, the temperature at which the enzymes were 

3 0 50% inactivated after 10 minutes (Tm) was calculated by 
extrapolation from the thermal inactivation plots and 
transition state theory was used to calculate 
activation energies for thermal inactivation (AG 1 ) . 
Table 7 shows the changes in AG* (AAG* ) and Tm for the 

35 combined mutant GAs relative to wild- type GA. These 

data clearly demonstrate that combining the individual 
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stabilizing mutations can cumulatively stabilize the 
enzyme . 

Saccharification analysis 

Figure 6 shows the results of saccharification 
5 analysis at 55° and 65°C for wild- type, S3 0P/G137A and 
S-S/S30P/G137A GAs using the industrial DE 10 
maltodextrin substrate Maltrin M100 (28% w/v) from 
Grain Processing Corporation. Complete conversion of 
28% w/v DE 10 maltodextrin to glucose would result in a 

10 1.71 M glucose syrup however, previous saccharification 
analyses in our laboratory have demonstrated that wild- 
type GA results in approximately 90% theoretical 
maximum glucose yield at 55°C (not shown) . At 55 6 C no 
significant difference in glucose production was 

15 observed between the wild-type and mutant enzymes. 
However, at 65°C the mutant GAs produced 8-10% more 
glucose than wild-type although none of the enzymes 
tested produced as much glucose as at 55°C probably due 
to thermal inactivation at the elevated reaction 

20 temperature. 

In summary these data show that the S30P/G137A 
double mutant enzyme was more stable than either single 
mutant GA when analyzed for resistance to irreversible 
thermal inactivation between 65°C and 80°C. The S- 

25 S/S30P combined GA mutant was also more stable than 
either the S30P or the S-S mutant GAs. The S- 
S/S30P/G137A combined mutant was the most stable GA 
variant constructed, particularly at temperatures above 
70°C when inactivated in a buffer system lacking mono- 

30 or polysaccharides. Saccharification analysis showed 
that the mutant enzymes performed better at elevated 
temperatures than wild-type GA. Importantly, none of 
the combined mutant GAs showed decreased enzyme 
activity when analyzed at 50°C. 

35 Discussion 

Sites of mutation 
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As described in Example 2, the mutations Asn20-*Cys 
and Ala27-»Cys form a disulfide bond between the C- 
terminus of a-helix one and an extended loop between ar- 
helices one and two. S3 OP and G137A were designed to 
5 stabilize the enzyme by reducing its conformational 

entropy of unfolding and are the most stabilizing in a 
series of proline substitution (Xaa->Pro) and Gly-»Ala 
mutations respectively. Ser30 is located at the second 
position of a type II 0-turn on an extended loop 

10 between a-helices one and two and Glyl37 is located in 
the middle of the fourth or-helix. 

It is of particular importance to note the 
positions of the S30P and the disulfide bond forming 
mutations. The disulfide bond is formed between 

15 positions 20 and 27; relatively close to position 30, 

The fact that both the disulfide bond forming mutations 
and the S3 OP stabilize GA suggests that this region of 
the enzyme is critical for irreversible thermal 
inactivation and may represent a region of local 

20 unfolding important for thermal inactivation. 

Additionally, previous investigators have suggested 
that a disulfide bond should not be engineered within 
four amino acids of a proline in primary sequence 
[Balaji et al, 1989] . This Example demonstrates that 

25 this rule is not absolute since thiol analysis showed 
that the disulfide bond was formed in the S-S/S30P and 
S-S/S30P/G137A combined mutants and thermal 
inactivation studies showed the stabilizing effects of 
the mutations were cumulative. 

30 Cumulative stabilization 

Previous work by Applicants has shown that 
combining two stabilizing mutations does not 
necessarily stabilize GA [Chen et al, 1996] . The 
present study, however, demonstrates that combining 

35 stabilizing mutations, even mutations very close to 

each other in the protein, can cumulatively stabilize 
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GA as measured by resistance to irreversible thermal 
inactivation. 

The S30P/G137A mutant showed more than additive 
stabilization at low temperatures (65°-70°C) , but less 
5 than additive stabilization at high temperatures 
(77.5°-80°C) (Figure 5A and Table 7) . At 80°C the 
inactivation rate for the S30P/G137A combined mutant 
was nearly identical to the S3 OP individual mutant 
protein. This indicates that both regions are very 

10 important for low temperature thermal inactivation, but 
at high temperatures inactivation became governed by 
other processes. 

It was somewhat surprising that combining the S3 OP 
with the disulfide bond forming mutations resulted in 

15 cumulative stabilization. This is not only because the 
engineered disulfide bond is so close to the engineered 
proline as is discussed above, but also because both 
are targeting the same region of the protein (ie: the 
extended loop between a-helices one and two) . It was 

20 expected that either the disulfide bond or S3 OP 
stabilized this region maximally, and further 
stabilization at this site would not result in a 
functionally more stable enzyme. As can be seen in 
Figure 5B, this was not the case. Combining the 

25 mutations resulted in roughly additive stabilization at 
all temperatures examined between 65° and 80 °C. 

The S-S/S30P/G137A combined mutant was no more 
stable than S30P/G137A GA at low temperatures (65°- 
70°C), but was slightly more stable at higher 

30 temperatures (75°-80°C) (Figure 5C and Table 7). 

Interestingly, the S-S/S30P GA is also more stable than 
S30P/G137A GA at high temperatures. Therefore, it 
appears that the introduced disulfide bond is 
particularly effective at stabilizing GA at high 

35 temperatures. 
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EXAMPLE 5 
INDUSTRIAL APPLICATION 

To determine whether the thermal stabilizing 
mutations: S30P/G137A and S-S/S30P/G137A would enhance 
5 GA performance under industrial conditions wild-type 
and mutant enzymes were subjected to high temperature 
saccharif ications (Figure 6) . Saccharif ication 
analysis showed that the mutant enzymes out -performed 
wild-type at 65°C but not at 55°C likely due to their 
10 increased stability. 
Conclusion 

The S30P/G137A double mutant cumulatively 
stabilized GA as demonstrated by decreased irreversible 
thermal inactivation rates relative to either 

15 individual mutant enzyme when analyzed between 65 °C and 
B0°C. Similarly, the S-S/S30P combined mutant also 
demonstrated cumulative stabilization. The S- 
S/S30P/G137A combined mutant was more stable than 
either of the "double" mutants, particularly at 

20 temperatures above 70°C. The S-S/S30P combined mutant 
had the same activity as wild-type and the S30P/G137A 
and S-S/S30P/G137A mutants increased enzyme activity by 
10-20% when assayed between 35° and 68 °C. The 
S30P/G137A and S-S/S30P/G137A mutant GAs decreased 

25 thermal inactivation rates approximately three fold 

relative to wild-type when inactivated in the presence 
of 1.71M glucose at 65°C. Additionally, at 55°C no 
difference in glucose yield was observed between these 
mutant GAs and wild-type for the saccharif ication of 

30 the industrial substrate Maltrin M100, whereas at 65°C 
the S30P/G137A and S-S/S30P/G137A GAs produced 8-10% 
more glucose than wild-type. 



35 



EXAMPLE 6 

MUTATIONS WITH INCREASED SELECTIVITY 

Interactions between substrates and charged 
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residues at subsites 1 and 2 of GA play a very 
important role in substrate specificity since the 
catalytic site is located between these sites. 
Therefore mutations were designed and analyzed to 
5 determine residues within these regions where mutations 
would increase selectivity of the enzyme reaction. In 
addition, several mutations that were designed to have 
thermostability were also screened for selectivity, as 
well as mutations designed to increase the pH optimum. 

10 Site-directed mutagenesis: Site-directed 

mutagenesis was performed as described herein above. 
The following mutagenic oligonucleotide primers were 
synthesized at the Iowa State University Nucleic Acid 
Facility: 5' -GGT CTC GGT GAG CCC AGG TTC AAT GTC GAT- 3 ' 

15 (Lysl08^Arg; SEQ ID No: 10), 5' -GGT CTC GGT GAG CCC ATG 
TTC AAT GTC GAT- 3 ' (Lysl08-*Met ; . SEQ ID No:ll), .5' -GAG 
GAC ACG TAC TGG AAC GGC AAC CCG-3 r (Tyr312^Trp; SEQ ID 
No: 12), and 5' -TAC CCT GAG GAC ACG TAC AAC GGC AAC GGC 
AAC TCG CAG GGC AAC CCG TGG TTC CTG TGC-3' (311-314 

20 Loop; SEQ ID No: 13), the underlined letters indicating 
the changed or added nucleotides. 
Results 

Enzyme kinetics 

As shown in table 11, the kinetic parameters Jc cat 

25 and for the hydrolysis of G 2 to G 7 as well as iG 2 in 
0.05 M acetate buffer, pH 4.4, at 45°C are given in 
Table 8. The 311-314 Loop mutant had k CAt values 50-80% 
for all a- (1,4) -linked substrates and only 30% for iG 2/ 
1C M values 50-75% for all substrates. The /c cat values for 

30 Glyl37^Ala/Ser30-»Pro GA are 10-30% more, generally, 
than that of wild-type GA for all substrates. The 
values of Glyl37-*Ala/Ser30-*Pro GA are about half to 
twofold for all the a- (1, 4) -linked substrates and 
essentially reached the wildtype level for iG 2 . The fc cat 

35 values for the GA engineered to carry the triple 

mutation, S-S/Glyl37^Ala/Ser30-»Pro, ranged from 80 to 
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120% generally for all substrates, and the values are 
30-80% for all substrates compared to wild-type GA. 
The Jc cat values for S-S GA are 85-110% for all 
substrates, and the S-S GA values are generally 90- 
5 110% for all substrates. However, the S-S GA iC M values 
are 140% for G 5 and 190% for G 6 . Values of Jc cat /J^ are 
75-105%, 60-110%, 60-110%, and 60-120% for the 
Tyr312-*Trp mutation, the combined Ser30->Pro/Glyl37-*Ala 
double mutation, the combined S-S/Ser30-*Pro/Glyl37-*Ala 
10 triple mutation, and the S-S engineered GA, 

respectively. The catalytic efficiencies for the 311- 
314 Loop GA are 85-120% for all the ot- (1, 4) -linked 
substrates, and only 50% for iG 2 , compared to wild- type 
GA. 

15 Table 8 shows the ratios of the catalytic 

efficiencies for G 2 to iG 2 for wild-type and mutant GAs. 
GAs engineered with the 311-314Loop mutation and 
Lysl08-*Arg mutation have the highest (240%) and the 
lowest (20%) catalytic efficiencies for or- (1,4)- over 

20 a- (1, 6) -linked substrates, respectively. The GAs 

engineered with the Tyr312-*Trp and S-S mutations show 
50% and 20% increases for this ratio, respectively. 
All other mutants had lower ratios, indicating poorer 
a- (1,4) -hydrolytic ability relative to a-(l,6)-hy- 

25 drolytic ability than wild-type GA. 
Mai tool igosaccharide hydrolysis 

GA engineered with the 311-314Loop mutation or 
with the S-S mutation had the highest average glucose 
yields (Figure 7) . The 311-314Loop GA had the lowest 

30 initial rates for glucose production (64%, 61%, and 82% 
compared to wild- type GA at 35, 45, and 55°C, 
respectively) due to a specific activity only 60% that 
of wild- type GA (data not shown) . Glucose 
concentrations decreased after reaching maximal values 

35 because of conversion to oligosaccharides. 
Glucose condensation reactions 
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IG 2 concentration profiles in 30% (w/v) glucose 
condensation reactions at 35, 45 , and 55°C were 
analyzed. GAs engineered with the Lysl08Arg mutation 
had the* highest and the 311-314Loop mutation as well as 
5 the S-S mutation the lowest equilibrium iG 2 

concentrations at all three temperatures. Tyr312->Trp, 
Ser30-»Pro/Glyl37-*Ala, and S-S/Ser30-*Pro/Glyl37->Ala GAs 
exhibited essentially the same equilibrium iG 2 
concentrations as wild-type GA. For all the other 
10 engineered thermostable GAs tested, Ser43 6-*Pro, S- 
S/Ser436-*Pro, S-S/Glyl37-»Ala, and 

Glyl37-*Ala/Ser436-*Pro, all reached higher equilibrium 
iG 2 concentrations than did wild- type GA. Table 9 shows 
the initial rates of iG 2 formation in 30% (w/v) glucose 

15 condensation reactions. S-S and 311-314Loop mutant GAs 
have the lowest initial rates at all three reaction 
temperatures tested. Lysl08->Arg mutant GA showed the 
highest initial rates among all the mutant GAs tested 
at all three reaction temperatures. All the tested 

20 thermostable GAs except Ser30->Pro/Glyl37-*Ala and S- 
S/Ser30-*Pro/Glyl37-*Ala had much higher initial rates 
than wild-type GA at 35°C, but they dropped to slightly 
higher or almost the same rate as wild- type GA at 55°C. 
The specificity for a- (1, 6) -linkage synthesis over a- 

25 (1,4) -linkage hydrolysis 

The ratio of the initial rate of iG 2 production in 
a 30% (w/v) glucose condensation reaction to that of 
glucose formation in 30% DE 10 maltodextrin hydrolysis 
was calculated to estimate the selectivity for the 

30 synthesis of ot- (1, 6) -linked products over the 

hydrolysis of ot- (1,6) -linked substrates. These 
iG 2 /glucose ratios and their relative ratios for wild- 
type and mutant GAs are given in Table 9. K108R and S- 
S mutants showed the highest and the lowest relative 

35 ratios among wildtype and all the mutant GAs at all 

reaction temperatures, respectively. Therefore, K108R 
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had more specificity for or- (1, 6) -linkages than a- (1,4)- 
linkages and S-S GA had more affinity for a-(l,4)- 
linkages than a- (1, 6) -linkages . The 311-314Loop GA 
also showed very low relative ratios at these three 
5 temperatures . 

EXAMPLE 7 

ADDITIONAL SELECTIVITY MUTATION ANALYSIS 

Utilizing the methods as set forth herein above, 

10 additional mutations were screened for selectivity as 
shown in Table 10 and Figures 8 and 9. 
Enzyme Kinetics: The kinetic parameters are seen in 
{ic cac and KJ for the hydrolysis of cr-1, 6-linked 
isomaltose and <*-l, 4 -linked maltooligodextrins (DP2-7) 

15 at 45°C and pH 4.4 are given in Table 10. Mutant Y175F 
was active. The k cat and K m values were 83-141% and 106- 
171%, respectively, that of wildtype for the different 
substrates tested and catalytic efficiencies were 69- 
102% that of wildtype. Mutant R241K was also active. 

20 Mutant S411G was highly active. The k cat and K m values 
were 93-129% and 83-203%, respectively, that of 
wildtype for the different substrates tested and 
catalytic efficiencies were 55-122% that of wildtype. 
Mutant S411A had a similar catalytic efficiency ratio 

25 as wildtype. Mutants Y116W, R241K, and S411G had 

decreased catalytic efficiency ratios compared to that 
of wildtype GA. 

DE 10 Maltodextrin Hydrolysis: At 55°C, the highest 
glucose yield was about 95% reached by engineered GA 

30 with mutant S411A at 216 hours compared to the wildtype 
yield of about 90% (Fig. 9) . All of the GAs, except 
S411A, reached their highest glucose yields rapidly. 
The glucose yield of S411A slowly increased for an 
extended period of time. The initial rates of glucose 

35 production at 55°C were 5 to 8 times higher that those 
at 35°C. 
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Glucose Condensation Reaction: Glucose condensation 
reactions were used to study the ability of wildtype 
and mutant GAs to synthesize isomaltose at high glucose 
concentrations (Figure 8) . The same concentrations of 
5 glucoamylases (2.64 /iM) were used as in the hydrolysis 
of DE 10 Maltodextrin. 

At 55°C, in spite of the different initial rates 
of isomaltose production for wildtype, R2 4 IK and Y175F, 
isomaltose production reached almost the same 

10 concentration at the last time point for these three 
mutant GAs (Figure 8) , indicating that the isomaltose 
production was close to equilibrium status. Isomaltose 
production for S411A and S411G was much lower than 
wildtype and almost linear as it was also at 35°C. 

15 Unexpectedly, isomaltose production for Y116W had a 
different (lower) equilibrium status compared to 
wildtype. The initial rates of isomaltose production 
at 55 °C were 5 to 7 times greater than those at 35 °C. 
R241K had a decreased initial rate of isomaltose 

20 production at 55 °C compared to that of wildtype, and it 
also had a lower increase (about 5 times) in the 
initial rate of isomaltose production from 35 °C to 
55 °C, compared to the wildtype increase (about 7 
times) . Y116W, Y175F, S411A and S411G had increased 

25 initial rates of isomaltose production or about 7, 6, 
and 5 times, respectively from 35°C to 55°C. 
Selectivity: The ratio of the initial rate of 
isomaltose production (from glucose condensation 
reactions) to that of glucose production (from 

30 hydrolysis of DE 10 maltodextrin) was calculated to 

evaluate selectivity for the synthesis of a- 1,6 -linked 
products versus the hydrolysis of a-1,4 linked 
substrates. This ratio represents the ability of a GA 
to synthesize isomaltose at a normalized level of DE 10 

35 maltodextrin hydrolytic activity. 

Mutants Y175F, S411A and S411G had a decreased 
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ratio of the initial rate of isomaltose production to 
that of glucose production to that of glucose 
production by 12%, 35% and 56% at 35°C, respectively, 
and a decreased ratio by 24%, 60% and 62% at 55°C, 
5 respectively, compared to wildtype. R241K had a very 
similar ratio to that of wildtype at both 35°C and 
55°C. 

EXAMPLE 8 

10 MUTATIONS TO PROVIDE pH OPTIMIZATION 

Utilizing the methods as set forth herein above, 
additional mutations, S411G, S411A, S411C, S411H, S411D 
were screened for increased pH optimum as shown in 
Figure 10 and Tables 11 and 12. 

15 Enzyme kinetics 

The kinetic parameters, k cat and K m , for the 
hydrolysis of a- 1, 4 -linked maltose and maltoheptaose 
and or-1, 6-linked isomaltose at 45°C and pH 4.4 are 
given in Table 11. Mutant S411G glucoamylase was 

20 highly active compared to wild- type, with an increased 
k cat and K m 13 - 30% and 11 - 59%, respectively, on the 
substrates tested. The catalytic efficiencies (k cac /KJ 
were 71 - 116% that of wild-type. Mutant S411A 
maintained 65 - 74% of wild-type catalytic efficiency 

25 with a slightly decreased k cat and a slightly increased 
K m . Mutant S411C maintained 54 - 73% of wild-type 
catalytic efficiency with a decrease in both the k cat • 
and K n values. Since mutant S411H and S411D had only 
about 6 - 12% of wild-type catalytic efficiency 

30 resulting from a seriously decreased k cat and an 

increased K n , the kinetic parameters for the hydrolysis 
of isomaltose were not determined. Only mutant S411H 
and S411D had large increases (5.5 to 7.5 kJ/mol) in 
the transition-state energy, A (AG), for the hydrolysis 

35 of maltose and maltoheptaose. The large increases of 

transition- state energy indicated that the introduction 
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of histidine or aspartic acid into position 411 
substantially destabilized the binding between GA and 
substrate in the transit ion- state. 
pH dependence of GA activity 
5 The kinetic parameters, Kat/Km and k cat , of the 

hydrolysis of maltose by wild-type and mutant 
glucoamylases at different pH values were calculated 
from initial rates obtained at low (smaller than 0.2 K m ) 
and high (higher than 10 K m ) concentrations of maltose. 

10 The effects of pH on the k cat /K m and k cat of maltose 

hydrolysis were used to determine the pK values (Table 
12) of both the free enzymes and the enzyme -substrate 
complexes. Although wild-type GA had a higher 
catalytic efficiency (k^/Kj than all of the mutant 

15 glucoamylases at all of the pH values tested, mutants 
S411G and S411A had higher k cat values than that of 
wild- type at some pH values. The uncomplexed and 
maltose -complexed S411H and S411D showed more narrow 
bell-shaped curves than that of wild-type. 

20 The effects of pH on the hydrolysis of 

maltoheptaose by wild- type, S411G and S411A GAs were 
measured to further investigate the change of pK values 
and optimum pH of enzyme- substrate complexes using a 
long- length substrate. Surprisingly, not only S411G, 

25 but also S411A were highly active compared to wild-type 
at the optimum pH. Wild-type GA pKl values (ionization 
of the catalytic base) were 2.77, 2.11, and 2.6 for the 
free enzyme, the maltose -complexed form, and the 
maltoheptaose -complexed form, respectively. The pK 2 

30 values (ionization of the catalytic acid) of wild-type 
were 5.80, 5.85, and 6.78 for the free enzyme, the 
maltose -complexed form, and the maltoheptaose -complexed 
form, respectively [Bakir et al . , 1993, Hiromi et al., 
1966, Sierks and Svensson, 1994]. Compared to wild- 

35 type, the S411G mutation increased the pKl of both the 
maltose -complexed form and the maltoheptaose-complexed 
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form by approximately 0.6 units, whereas S411G had no 
effect on the pK2 of either enzyme -substrate complexes 
and only had a minor effect on the pKl and pK2 of the 
free enzyme. The combined effect of S411G on pKl and 
5 pK2 was an increased optimum pH of both the maltose - 
complexed form and the maltoheptaose-complexed form by 
approximately 0.3 units. 

The S411G mutation, however, had no effect on the 
optimum pH of the free enzyme. S411A and S411C had very 

10 similar effects on the pH dependence of maltose 

hydrolysis. S411A and S411C increased the pK x of the 
free enzyme and the maltose -complexed forms by 0.3 - 
0.5 and 1.21 units, respectively. Surprisingly, S411A 
and S411C also increased the pK 2 of the maltose- 

15 complexed form by approximately 0.5 units. In 
addition, S411A increased the pKj and pK 2 of the 
maltoheptaose-complexed form by 1.31 and 0.4 units, 
respectively. S411H increased the pK x of the free 
enzyme and maltose -complexed form by 0.33 and 1.47 

20 units, respectively; however, it decreased the pK 2 of 
the free enzyme and the maltose -complexed form by 0.79 
and 1.16 units, respectively. S411D increased the pKl 
of the free enzyme and the maltose -complexed form by 
0.36 and 1.23 units, respectively. S411D also 

25 decreased the pK 2 of the maltose -complexed form by 0.32 
units. For wild- type, S411G, and S411A GAs, the values 
of pK 1# pK 2/ and pH opt for the maltoheptaose-complexed 
forms were higher than those of the corresponding 
maltose -complexed forms by approximately 0.5, 0.9 and 

30 0.7 units, respectively. For S411G and S411A, the 

increases in pH optimum (compared to that of wild- type) 
obtained using the long- length substrate 
(maltoheptaose) were almost the same as that obtained 
using the short -length substrate (maltose) . 

35 All five mutants at position 411 showed a shift of 

0.15 to 0.87 units in the optimum pH of the enzyme- 
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substrate complex compared to wild- type (Table 12) , 
mainly due to increased pK x values. Compared with other 
mutants, S411A was the best performing pH mutant. 
S411A increased the optimum pH by 0.84 units while also 
5 maintaining a high level of both catalytic activity 
(k cat ) and catalytic efficiency (k cat /K m ) . 
The hydrolysis of maltodextrin 10 

The hydrolysis of 28% (w/v) maltodextrin was used 
to study the pH dependence of GA activity at a high 

10 concentration of a long-length substrate. Maltodextrin 
10 is a mixture of maltodextrin with an average (and 
major) degree of polymerization of 10. The production 
of glucose by wild-type and S411A glucoamylases during 
the hydrolysis of maltodextrin 10 at 11 different pH 

15 values was determined, and used to calculate the 

initial rates of glucose production at different pH 
values (Figure 10) . The production of glucose 
increased following a hyperbolic curve. S411A had 
higher initial rates of glucose production than wild- 

20 type when the pH values were above 6.6 (Figure 10) . 

Throughout this application, various 
publications are referenced by author and year and 
patents listed by number. Full citations for the 
publications are listed below. The disclosures of 

25 these publications and patents in their entireties are 
hereby incorporated by reference into this application 
in order to more fully describe the state of the art to 
which this invention pertains . 

The invention has been described in an 

30 illustrative manner, and it is to be understood that 
the terminology which has been used is intended to be 
in the nature of words of description rather than of 
limitation. 

Obviously, many modifications and variations of 
35 the present, invention are possible in light of the 

above teachings. It is, therefore, to be understood 
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that within the scope of the appended claims, the 
invention may be practiced otherwise than as 
specifically described. 
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TABLE 1 



Changes in AG* and Tm for the mutant GAs relative to wild-type. 





AAG* 


ATm 


GA form 


(kJ/mol.) 


(°C) 


Ser30->Pro 


1.6 


1.7 


Asp345-*Pro 


0.5 


0.4 


Glu408->Pro 


-7.2 


-6.7 



TABLE 2 



Summary of DTNB-titratable sullhydryl groups in wild-type and mutant GA 
with or without DTT reducing 



[SH]/molecule No. of disulfide bonds* 

DTT+ DTT- 

WT 8 . 6 
A27C/N20C 10.9 
A471C/T72C 10.4 

* No. of disulfide bonds= ([SH]/molecule (DTT+)-[SH]/molecule (DTT-))/2 



Enzymes 



0.9 4 
0.9 5 
1.3 5 
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TABLE 3 



Catalytic properties of wild type and mutant GAs. 



GAform Specific Activity K m kcjfl ffcai^m 

(IU/mgGA) (mM) (s ) (s mM ) 

— m vat 



WT* 


20.6±0 2 


0.72±0.03 


8.67±0.17 


12.0 


A27C 


14.9*1.1 


0.86±0.1t 


8.02±0.45 


9.3 


N20C 


8.1 ±0.5 


0.70±0.05 


3.97±0.12 


5.7 


A27C/N20C 


18.3±0.7 


0.90±0.08 


9.61±0.40 


J0.7 


A471C/T72C 


22.7±1.5 


0.87±0.07 


I0.17±0.40 


11.6 


A27C/N20C 
/S436 


22.5±1.8 


n/d c 


N/D 


N/D 


A27C/N20C 
/G137A 


24.2±0.8 


N/D 


N/D 


N/D 


G137A/S436P 


25.0±0.9 


N/D 


N/D 


N/D 



* Produced in shaking flasks 

standard error 
c not determined 
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TABLE 4 



Activation parameters for irreversible thermoinactivation of wild-type (WT) 
and mutant G As at pH 4,5. 



GAform AH* AS* AG * (65°C) AG*(75°C) 

(kJ/mol) (J/molK) (kJ/mol) (kJ/mol) 



WT* 


366± I b 


769± 4 


105.7 


98.0 


A27C 


370±15 


780±44 


106.3 


98.5 


N20C 


324±1 1 


654*33 


102.7 


96.2 


A27C/N20C 


342±16 


694±46 


107.2 


100.2 


A471C/T72C 


365± 9 


768±26 


106.0 


98.3 


A27C/N20C 
/S436P 


352± 6 


724±18 


107.9 


100.7 


A27C/N20C 
/G137A 


362± 1 


751± 2 


108.4 


100.9 


G137A/S436P 


362*20 


752±57 


107.7 


100.2 


S436P c 


351± 8 


723±24 


106.2 


99.0 


G137A<* 


330± 6 


66I±I7 


106.5 


99.9 



* produced by shaking flask 

standard error 
c UetaL, 1996 
d Chene/a/., 1996 
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Table 5. Specific activities of wild-type and mutant GAs 





Specific activity* 


GA form 


(IU/mg) 


Wild-type 


21.1±0.1 


S30P/GIyl37A 


24.0±1.2 


S-S/S30P 


2l.2±0.5 


S-S/S30P/G137A 


24.5 ±0.2 



* Standard deviation resulting from three or more assays 



Table 6. Thiol analysis of wild-type and mutant GAs 



GA form 


[Protein] (nM) 


[SH] (nM) a 


[SH]/[Proteinl* 


Wild-type 


10 


8 


0.8 


S-S/S30P 


10 


11 


1.1 


S-S/S30P/G137A 


10 


13 


1.3 



a Average of duplicate analysis 



Table 7. Changes in free energies for thermal inactivation (AAG*) and 
temperatures at which the enzyme is 50% inactivated after 10 minutes (ATm) 
relative to wild-type GA 





AAG** 




GA form 


(kJ/mol) 


ATm(°C) 


8301* 


1.6 


1.7 


G137A C 


0.8 


1.2 


S-S d 


1.2 


1.4 


S30P/G137A 


4.5 


3.5 


S-S/S30P 


3.5 


3.2 


S-S/S30P/G137A 


4.4 


3.9 



1 Calculated at 65°C 
b From Allen et. a/. 8 
c From Chen et. aL 6 
d From Li et. al? 
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TABLE 9 

Initial rates of glucose and isomaltose productions in the hydrolysis 
of 30% (w/v) maltodextrin M100 and 30% (w/v) glucose condensations, 
respectively, and their relative ratios for wild- type and mutant glucoamylases 
at 35°C, 45°C, and 55°C. 



Initial rates Ratios 
Enzymes Glucose* (Gl) Isomaltose* (iG2> Ratios Relative 

(pg/mL-h)xlO~ 3 (pig/mL*h)x!0 3 (iG 2 /Gl)xl0 € ratios 



Wild-type 


21. 


5 


±0. 


6° 


289 


± 


5 


13. 


5 


1.00 


LyslOSArg 


22. 


0 


+0. 


4 


969 


± 


12 


44. 


1 


3.27 


Tyr312Trp 


17. 


9 


±0. 


5 


294 


± 


4 


16. 


4 


1.21 


311-314Loop 


13. 


8 


±0. 


4 


128 


± 


3 


9. 


3 


0.69 


Ser30Pro/Glyl37Ala 


27. 


7 


±0. 


4 


298 


± 


6 


10. 


8 


0.80 


S-S/Ser30Pro/Glyl37Ala 


30. 


1 


±0. 


5 


245 


± 


6 


8. 


2 


0.60 


S-S 


31. 


8 


±0. 


6 


135 


± 


3 


4. 


2 


0.31 


Ser436Pro 


29. 


9 


±0. 


6 


903 


± 


12 


30. 


2 


2.23 


S-S/Ser436Pro 


31. 


7 


±0. 


5 


824 


± 


12 


26. 


0 


1.92 


S-S/Glyl37Ala 


35. 


2 


±0. 


6 


982 


± 


15 


27. 


9 


2.07 


Glyl37Ala/Ser436Pro 


36. 


6 


+0. 


7 


776 


± 


10 


21. 


2 


1. 57 


45°C 






















Wild- type 


66. 


2 


±2. 


2 


3880 


± 


60 


58. 


7 


1.00 


Lysl08Arg 


50. 


.2 


±2. 


.0 


6420 


± 


110 


128 




2.18 


Tyr312Trp 


52. 


6 


±2. 


, 1 


3360 


± 


60 


63. 


9 


1.09 


311-314Loop 


40. 


.4 


±1. 


.8 


1430 


± 


40 


35. 


3 


0.60 


Ser30Pro/Glyl37Ala 


76. 


.3 


±2, 


, 7 


3690 


± 


70 


48. 


4 


0.83 


S- S /Ser3 0Pro/Glyl37Ala 


84. 


.3 


±3. 


,0 


3520 


± 


60 


41. 


7 


0.71 


S-S 


86.3 


±3, 


.3 


963 


± 


28 


11. 


.2 


0.19 


S5°C 






















Wild-type 


156 


±3 




4890 


± 


80 


31 


.3 


1.00 


Lysl08Arg 


101 


+ 1 




8200 


± 


120 


81 


.0 


2.59 


Tyr312Trp 


110 


±2 




4440 


+ 


70 


40 


.5 


1.29 
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TABLE 9 (continued) 



311-314Loop 


12B 


±2 


1890 


± 


50 


14. 


8 


0 


.47 


Ser30Pro/Glyl37Ala 


157 


±3 


7200 


± 


110 


«<■ 


8 


1. 


.47 


S-S/Ser30Pro/Glyl37Ala 


167 


±3 


5690 


± 


100 


34. 


1 


1. 


.09 


S-S 


164 


±3 


1230 


± 


40 


7. 


5 


0. 


.24 


Ser436Pro 


218 


±3 - 


4710 


± 


80 


21. 


6 


0. 


.69 


S-S/Ser436Pro 


ND d 




5130 


± 


100 


ND 




ND 




SS/Glyl37Ala 


225 


±3 


5720 


± 


100 


25. 


4 


0. 


81 


Glyl37Ala/Ser436Pro 


208 


±3 


ND 






ND 




ND 





Samples were taken from 30% (w/v) M100 hydrolysis reactions in 0.05 M NaOAc 
buffer, pH 4.4; glucose concentrations were determined by glucose oxidase 
method. 

b Samples were taken from 30% (w/v) glucose condensation reactions in 0.05 M 

NaOAc buffer, pH 4.4; isomaltose concentrations were determined by HPTLC. 

Standard error 
d 

Not determined 
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TABLE 13 



Increases in free energies for thermal inactivation (A AG*) relative 
to wild-type GA calculated at 65° C 





AAG* 


GA form 


(kJ/mol) 


S436P 


0.5 


S30P 


1.6 


G137A 


0.8 


S-S 


1.2 


S-S/S436P 


2.2 


G137A/S436P 


2.0 


S30P/G137A 


4.5 


S-S/S30P 


3.5 


S-S/G137A 


2.7 


S-S/S30P/G137A 


4.4 



* A AG* greater than zero indiactcs increased thermostability 



TABLE 14 

Decrease in the relative ratio of initial rate of isomaltose formation from 30% (w/v) glucose 
condensation reactions to that of glucose formation in 30% (w/v) maltodextrin Ml 00 hydrolysis 
reactions. 



GAform 


Relative ratios" 


Wild-type 


1.00 


S-S 


0.24 


S30P 


0.77 


G137A 


0.54 


Y175F 


0.76 


300Loop 


0.47 


S411A 


0.40 


S411G 


0.38 


S436P 


0.70 


S-S/G137A 


0.81 


G121A/S411G 


0.44 



All the above reactions were carried out in 0.05M sodium acetate buffer, pH4.4, at 55° C. 
'Ratios lower than 1.00 indicate increased specificity for a-(l,4) over a-(l,6)-linked substrates. 
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TABLE 15 



Increase in the optimum pH of the enzyme-substrate complex of mutant glucoamylases for 
hydrolysis of maltose at 45° C compared to that of wild-type. 



G A form pH^ Increase" 



S411G 0.26 

S411A 0.84 

S411C 0.86 

S411H 0.15 

S411D 0.46 



*The pH optimum of the enzyme-substrate complex of wildtype glucoamylase for hydrolysis of 
maltose at 45° C was pH 3.98. 



WO 98/03639 



-69- 



PCTAJS97/12983 



REFERENCES 

Ahearn and Klibanov. The Mechanism of Irreversible 
Enzyme Inactivation at 100°C. Science, 228, 1280 
(1985) . 

Ahearn, et al . Control of Oligomeric Enzyme 
Thermostability by Protein Engineering. Proc. Nat. 
Acad. Sci. U. S. A., 84, 675 (1987). 

Aleshin, et al., 1992. Crystal structure of 
glucoamylase from Aspergillus awamori var. X100 to 2.2- 
A resolution. J. Biol. Chem. 267:19291-19298. 

Aleshin, et al., 1994. Refined structure for the 
complex of acarbose with glucoamylase from Aspergillus 
awamori var. X100 to 2. 4 -A resolution. J. Biol. Chem. 
269 :15631-15639. 

Aleshin et al., 1994, J. Mol . Biol, 238:575-591 

Aleshin, et al., 1996. Crystallographic complexes of 
glucoamylase with maltooligosaccharide analogs: 
relationship of stereochemical distortions at the 
nonreducing end to the catalytic mechanism. 
Biochemistry. 35: 8319-8328. 

Argos, et al . Thermal Stability and Protein Structure. 
Biochemistry, 18, 5698 (1979) . 

Bakir et al (1993) Protein Eng. 6:939-946 

Balaji, et al 1989. Modification of protein stability 
by introduction of disulfide bridges and prolines; 
Geometric criteria for mutation sites. Biochem. 
Biophys. Res. Commun. 160:109-114. . 

Banks, W. and Greenwood, C. J. 1971. The specific 
estimation of glucose using glucose oxidase. StMrke. 
23:222-228. 

Boel, et al. Glucoamylases Gl and G2 from Aspergillus 
niger Are Synthesized from Two Different but Closely 
Related mRNAs. EMBO J. , 3, 1097 (1984). 

Bradford. A Rapid and Sensitive Method for the 
Quantitation of Microgram Quantities of Protein 
Utilizing the Principle of Protein— Dye Binding. Anal. 
Biochem., 72, 248 (1976). 

Chen, et al (1996) . Effect of replacing helical 
glycine residues with alanine on reversible and 
irreversible stability and production of Aspergillus 
awamori glucoamylase. Protein Eng. 9:499-505. 



WO 98/03639 



PCT/US97/12983 



-70- 



Chen, et al (1995) Protein Eng., 8, 575-582. 

Chen, et al (1994a) . Increased thermostability of 
Asnl82-»Ala mutant Aspergillus awamori glucoamylase . 
Biotechnology and Bioengineering . 43:101-105. 

Chen, et al (1994b) Biochem. J., 301, 275-281. 

Clarke and Svensson. Identification of an Essential 
Tryptophanyl Residue in the Primary Structure of 
Glucoamylase G2 from Aspergillus niger. Carlsberg Res. 
Commun., 49, 559 (1984). 

Clarke and Fersht, (1993) Biochemistry, 32, 4322-4329. 

Cole, et al. Stable Expression of Aspergillus awamoi 
Glucoamylase in Distiller's Yeast. Bio/ Technol., 6, 
All (1988) . 

Coutinho and Reilly, (1994a) Protein Eng., 7, 393-400. 

Coutinho and Reilly, 1994b. Structural similarities in 
glucoamylases by hydrophobic cluster analysis. Protein 
Eng. 7:749-760. 

Fierobe et al., 1996. Mutational modulation of 
substrate bond-type specificity and thermostability of 
glucoamylase from Aspergillus awamori by replacement 
with short homologue active site sequences and 
thiol/disulfide engineering. Biochemistry. 35:8696- 
8704. 

Frandson et al . , 1994. Biochemistry 33:13808-13816. 

Garfin, (1990) In: Guide to protein purification. M. P. 
Deutscher, Ed. Methods in Enzymology. vol. 182. 
Academic Press, San Diego, CA. pp. 425-441. 

Gunnarsson, et al. Structural Studies on the O-Glycos- 
idically Linked Carbohydrate Chains of Glucoamylase Gl 
from Aspergillus niger. Eur. J. Biochem. , 145, 463 
(1984) . 

Habeeb, (1972). Methods in Enzymology. vol. 25. 
Academic Press, San Diego, CA. pp. 457-464. 

Harris et al., 1993. Refined structure for the complex 
of l-deoxynojirimycin with glucoamylase from 
Aspergillus awamori var. X100 to 2.4-A resolution. 
Biochemistry. 32:1618-1626. 

Hiromi et al . , 1983, Mol. Cell. Biochem, 51:79-95 
Hiromi et al. (1966) J. Biochem 59:469-475 



WO 98/03639 



PCT/US97/12983 



-71- 



Hurley, et al . , 1992. Flexible -geometry 
conformational energy maps for the amino acid residue 
preceding a proline. Biopolymers. 32:1443-1446. 

Imanaka, et al. A New Way of Enhancing the 
Thermostability of Proteases. Nature, 324, 695 (1986). 

Innis, et al. Expression, Glycosylation, and Secretion 
of an Aspergillus Glucoamylase by Saccharomyces 
cerevisiae. Science, 228, 21 (1985) . 

Itoh, et al. Nucleotide Sequence of the Glucoamylase 
Gene GLU1 from Saccharomycopsis fibuligera. J. 
Bacteriol., 169, 4171 (1987). 

Klibanov. Stabilization of Enzymes Against Thermal 
Inactivation. Adv. Appl . Microbiol., 29, 1 (1983). 

Knowles, et al. Cellulase Families and Their Genes. 
Trends Biotechnol., 5, 255 (1987). 

Kunkel . Rapid and Efficient Site-Specific Mutagenesis 
Without Pheno type Selection. Proc. Nat. Acad. Sci. C7. 
S. A . , 82, 448 (1985) . 

Kunkel, et al 1987. Rapid and efficient site-specific 
mutagenesis without phenotypic selection. Methods 
Enzymol . 154 : 367-382 . 

Lee, et al . Pilot Plant Production of Glucose with 
Glucoamylase Immobilized to Porous Silica. Biotechnol. 
Bioeng., 16, 1507 (1976). 

Leatherbarrow and Fersht . Protein Engineering. 
Protein Eng., 1, 7 (1986). 

Matthews, et al., 1987. Enhanced protein 
thermostability from site-directed mutations that 
decrease the entropy of unfolding. Proc. Natl. Acad. 
Sci. U. S. A. 84:6663-6667. 

MacArthur and Thornton, 1991. Influence of proline 
residues on protein conformation. J. Mol. Biol. 
218:397-412. 

Masumura, et al (1989) Nature, 342, 291-293. 

Matsumura and Aiba. Screening for Thermostable Mutant 
of Kanamycin Nucleotidyltransferase by the Use of a 
Transformation System for a Thermophile, Bacillus 
stearothermophilus. J. Biol. Chem. , 260, 15298 (1985). 

Matsumura, et al. A Cumulative Effect of Intragenic 
Amino Acid Replacements on the Thermostability of a 



WO 98/03639 



PCI7US97/12983 



-72- 



Protein. Nature, 323, 356 (1986). 

Matsuura, et al. Structure and Possible Catalytic 
Residues of Taka-Amylase A. J. Biochem. , 95, 697 
(1984) . 

Mcllvane (1921) Biochem. J. 49:183-6. 

Meagher et al., 1989, Biotechnol Bioeng. 34:681-688 

Munch and Tritsch. Irreversible Thermoinactivation of 
Glucoamylase from Aspergillus niger and 
Thermostabilization by Chemical Modification of 
Carboxyl Groups. Biochim. Biophys. Acta, 1041, 111 
(1990) . 

Nikolov, et al . Kinetics, Equilibria, and Modeling of 
the Formation of Oligosaccharides from D-Glucose by 
Aspergillus niger Glucoamylases I and II. Biotechnol. 
Bioeng., 34, 694 (1989). 

Nosoh and Sekiguchi. Protein Engineering for 
Thermostability. Biocatal., 1, 257 (1988). 

Nunberg, et al. Molecular Cloning and Characterization 
of the Glucoamylase Gene of Aspergillus awamori. Mol. 
Cell. Biol., 4, 2306 (1984). 

Pakula and Sauer. Genetic Analysis of Protein 
Stability and Function. Arm. .Rev. Genet., 23, 289 
(1989) . 

Perry and Wetzel. Disulfide Bond Engineered into T4 
Lysozyme: Stabilization of the Protein Toward Thermal 
Inactivation. Science, 226, 555 (1984) . 

Pollitt and Zalkin, (1983). J. Bacteriology, 153, 27- 
32. 

Rabbo and Terkildsen (1960). Scandinav. J. & Lab. 
Investigation 12:402-407 

Ramachandran, et al . , 1963. Stereochemistry of 
polypeptide chain configurations. J". Mol . Biol. 7:95- 
99. 

Robyt & Mukerjea (1994) Carbohydr. Res. 251:187-202 

Savel'ev and Firsov. Carboxyl Groups at the Active 
Site of Glucoamylase from Aspergillus awamori. 
Biochemistry (U. 5. S. R.), 47, 1365 (1982)'. 

Savel'ev, et al. Study of the Active Center of 
Glucoamylase from Aspergillus awamori. Biochemistry 



WO 98/03639 



PCT/US97/12983 



-73- 



(U. S. S. R.) , 47, 330 (1982) . 

Semimaru, et al (1995) Appl. Environ. Microbiol. , 61, 
2885-2990. 

Sierks & Svensson (1994) Protein Eng. 7:1479-1480 

Sierks, et al . Site-Directed Mutagenesis at the Active 
Site Trpl20 of Aspergillus awamori Glucoamylase . 
Protein Eng., 2, 621 (1989). 

Sierks, et al . Catalytic Mechanism of Fungal Glucoamy- 
lase as Defined by Mutagenesis of Aspl76, Glul79, and 
Glu 180 in the Enzyme from Aspergillus awamori. 
Protein Eng., 3, 193 (1990). 

Smith, et al., 1985. Measurement of protein using 
bicinchoninoc acid. Anal Biochem. 150:76-85. 

Stoffer, et al., 1995. Refined structure for the 
complex of D- gl uco- dihydroacarbose with glucoamylase 
from Aspergillus awamori var. X100 to 2.2 A resolution: 
dual conformations for the extended inhibitors bound to 
the active site of glucoamylase. FEBS Letters. 358:57- 
61. 

Suzuki, et al . , 1987. A strong correlation between 
the increase in the number of proline residues and the 
rise in thermostability of five Bacillus oligo-1,6- 
glucosidases. Appl. Microbiol. Biotechnol . 26:546- 
551. 

Suzuki, et al . , 1991. A hyperthermostable pullulanase 
produced by an extreme thermophile, Bacillus 
flavocaldarius KP 1228, and evidence for the proline 
theory of increasing protein thermostability. Appl. 
Microbiol. Biothnol . 34:707-714. 

Suzuki, 1989. A general principle of increasing 
protein thermostability. Proc. Jpn. Acad. Ser B. 
65:146-148. 

Suzuki, et al . Site-Directed Mutagenesis Reveals 
Functional Contribution of Thr218, Lys220, and Asp304 
in Chymosin. Protein Eng., 4, 69 (1990). 

Svensson. Regional Distant Sequence Homology Between 
Amylases, a-Glucosidases and Transglycosylases . FEBS 
Lett., 230, 72 (1988). 

Svensson, et al . Characterization of Two Forms of 
Glucoamylase from Aspergillus niger. Carlsberg Res. 
Commun. , 47, 55 (1982) . 



WO 98/03639 



PCT/US97/12983 



-74- 



Svensson, et al. The Complete Amino Acid Sequence of 
the Glycoprotein, Glucoamylase Gl f from Aspergillus 
niger. Carlsberg Res. Commun. , 48, 529 (1983). 

Svensson, et al. Sequence Homology Between Putative 
Raw-Starch Binding Domains from Different Starch- 
Degrading Enzymes. Biochem. J, Lett., 264, 309 (1989). 

Svensson, et al. Identification of Carboxylic Acid 
Residues in Glucoamylase G2 from Aspergillus niger That 
Participate in Catalysis and Substrate Binding. Eur. 
J. Biochem., 188, 29 (1990). 

Tanaka et al., 1983, J. Biochem., 93:1037-1043 

Watanabe, et al., 1994. Multiple proline substitutions 
cumulatively thermostabilize Bacillus cereus ATCC7064 
oligo-l,6-glucosidase. Eur. J". Biochem. 226:277-283. 

Wasserman. Thermostable Enzyme Production. Food 
Technol., 38, 78 (1984). 

Wetzel, R. (1987) Trends Biochem. Sci., 12, 478-482. 
Williamson, et al (1992) Biochem. J". , 282, 423-428. 



WO 98/03639 



PCT/US97/12983 



-75- 



SEQUKNCE LISTING 



(1) GENERAL INFORMATION: 

(i) APPLICANT: Allen, Martin 

Fang, Tsuei-Yun 
Li , Yuxing 
Liu, Hsuan- Liang 
Chen, Hsui-Mei 
Coutinho , Pedro 
Hanzatko, Richard 
Ford, Clark 

(ii) TITLE OF INVENTION: PROTEIN ENGINEERING OF GLUCOAMYLASE TO 
INCREASE pH OPTIMUM, SUBSTRATE SPECIFICITY AND 
THERMOSTABILITY 

(iii) NUMBER OF SEQUENCES: 12 

(iv) CORRESPONDENCE ADDRESS: 

(A) ADDRESSEE: Kohn & Associates 

(B) STREET: 30500 Northwestern Hwy. 

(C) CITY: Farmington Hills 

(D) STATE: Michigan 

(E) COUNTRY: US 

(F) ZIP: 48334 

(v) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS/MS-DOS 

(D) SOFTWARE: Patentln Release #1.0, Version #1.30 

(vi) CURRENT APPLICATION DATA: 

(A) APPLICATION NUMBER: 

(B) FILING DATE: 

(C) CLASSIFICATION: 

(viii) ATTORNEY/ AGENT INFORMATION: 

(A) NAME: Kohn, Kenneth I. 

(B) REGISTRATION NUMBER: 30,955 

(C) REFERENCE / DOCKET NUMBER: 0812.00001 

(ix) TELECOMMUNICATION INFORMATION: 

(A) TELEPHONE: (248) 539-5050 

(B) TELEFAX: (248) 539-5055 



(2) INFORMATION FOR SEQ ID NO:l: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 616 amino acids 

(B) TYPE: amino acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: protein 

(iii) HYPOTHETICAL: NO 

(vi) ORIGINAL SOURCE: 

(A) ORGANISM: Aspergillus 
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(xi) SEQUENCE DESCRIPTION: SEQ ID NO : 1 : 

Ala Thr Leu Asp Ser Trp Leu Ser Aan Glu Ala Thr Val Ala Arg Thr 
1 5 10 15 

Ala He Leu Asn Asn He Gly Ala Asp Gly Ala Trp Val Ser Gly Ala 
20 25 30 

Asp Ser Gly He Val Val Ala Ser Pro Ser Thr Asp Asn Pro Asp Tyr 
35 40 45 

Phe Tyr Thr Trp Thr Arg Asp Ser Gly Leu Val Leu Lys Thr Leu Val 
50 55 60 

Asp Leu Phe Arg Asn Gly Asp Thr Ser Leu Leu Ser Thr He Glu Asn 
65 70 75 80 

Tyr He Ser Ala Gin Ala He Val Gin Gly He Ser Asn Pro Ser Gly 
85 90 95 

Asp Leu Ser Ser Gly Ala Gly Leu Gly Glu Pro Lys Phe Asn Val Asp 
100 105 110 

Glu Thr Ala Tyr Thr Gly Ser Trp Gly Arg Pro Gin Arg Asp Gly Pro 
115 120 125 

Ala Leu Arg Ala Thr Ala Met He Gly Phe Gly Gin Trp Leu Leu Asp 
130 135 140 

Asn Gly Tyr Thr Ser Thr Ala Thr Asp He Val Trp Pro Leu Val Arg 
145 150 155 160 

Asn Asp Leu Ser Tyr Val Ala Gin Tyr Trp Asn Gin Thr Gly Tyr Asp 
165 170 175 

Leu Trp Glu Glu Val Asn Gly Ser Ser Phe Phe Thr He Ala Val Gin 
180 185 190 

His Arg Ala Leu Val Glu Gly Ser Ala Phe Ala Thr Ala Val Gly Ser 
195 200 205 

Ser Cys Ser Trp Cys Asp Ser Gin Ala Pro Glu He Leu Cys Tyr Leu 
210 215 220 

Gin Ser Phe Trp Thr Gly Ser Phe He Leu Ala Asn Phe Asp Ser Ser 
225 230 235 240 

Arg Ser Gly Lys Asp Ala Asn Thr Leu Leu Gly Ser He His Thr Phe 
245 250 255 

Asp Pro Glu Ala Ala Cys Asp Asp Ser Thr Phe Gin Pro Cys Ser Pro 
260 " 265 270 

Arg Ala Leu Ala Asn His Lys Glu Val Val Asp Ser Phe Arg Ser He 
275 280 285 

Tyr Thr Leu Asn Asp Gly Leu Ser Asp Ser Glu Ala Val Ala Val Gly 
290 295 300 

Arg Tyr Pro Glu Asp Thr Tyr Tyr Asn Gly Asn Pro Trp Phe Leu Cys 
305 310 315 320 

Thr Leu Ala Ala Ala Glu Gin Leu Tyr Asp Ala Leu Tyr Gin Trp Asp 
325 330 335 
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Lys Gin Gly Ser Leu Glu Val Thr Asp Val Ser Leu Asp Phe Phe Lys 
340 345 350 

Ala Leu Tyr Ser Asp Ala Ala Thr Gly Thr Tyr Ser Ser Ser Ser Ser 
355 360 365 

Thr Tyr Ser Ser lie Val Asp Ala Val Lys Thr Phe Ala Asp Gly Phe 
370 375 380 

Val Ser He Val Glu Thr His Ala Ala Ser Asn Gly Ser Met Ser Glu 
385 390 395 400 

Gin Tyr Asp Lys Ser Asp Gly Glu Gin Leu Ser Ala Arg Asp Leu Thr 
405 410 * 415 

Trp Ser Tyr Ala Ala Leu Leu Thr Ala Asn Asn Arg Arg Asn Ser Val 
420 425 430 

Val Pro Ala Ser Trp Gly Glu Thr Ser Ala Ser Ser Val Pro Gly Thr 
435 440 445 

Cys Ala Ala Thr Ser Ala He Gly Thr Tyr Ser Ser Val Thr Val Thr 
450 455 460 

Ser Trp Pro Ser He Val Ala Thr Gly Gly Thr Thr Thr Thr Ala Thr 
465 470 475 , 480 

Pro Thr Gly Ser Gly Ser Val Thr Ser Thr Ser Lys Thr Thr Ala Thr 
485 490 495 

Ala Ser Lys Thr Ser Thr Ser Thr Ser Ser Thr Ser Cys Thr Thr Pro 
500 505 ~ 510 

Thr Ala Val Ala Val Thr Phe Asp Leu Thr Ala Thr Thr Thr Tyr Gly 
515 520 525 

Glu Asn He Tyr Leu Val Gly Ser lie Ser Gin Leu Gly Asp Trp Glu 
530 535 540 

Thr Ser Asp Gly lie Ala Leu Ser Ala Asp Lys Tyr Thr Ser Ser Asp 
545 550 555 560 

Pro Leu Trp Tyr Val Thr Val Thr Leu Pro Ala Gly Glu Ser Phe Glu 
565 570 575 

Tyr Lys Phe He Arg He Glu Ser Asp Asp Ser Val Glu Trp Glu Ser 
580 585 590 

Asp Pro Asn Arg Glu Tyr Thr Val Pro Gin Ala Cys Gly Thr Ser Thr 
595 600 605 

Ala Thr Val Thr Asp Thr Trp Arg 
610 615 

(2) INFORMATION FOR SEQ ID NO: 2: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 7 amino acids 

(B) TYPE: amino acid 

(C) STRANDEDNESS : single 
<D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: peptide 
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(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2 

Asn Gly Asn Gly Asn Ser Gin 
1 5 

(2) INFORMATION FOR SEQ ID NO: 3: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 31 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE : other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3 
CAGAGTCCGC GCCCGGCACC CAAGCACCGT C 
(2) INFORMATION FOR SEQ ID NO:4: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 33 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = B PRIMER 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4 
AAGTCCAGCG ACACAGGTGT GACCTCCAAC GAC 
(2) INFORMATION FOR SEQ ID NO: 5: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 32 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO : 5 
CGAGCGGAAA GCTGCGGGCC ATCAGACTTG TC 
(2) INFORMATION FOR SEQ ID NO: 6: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 32 base pairs 

(B) TYPE: nucleic acid 
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(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER" 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:6: 
CGTACTGCCA TCCTGTGTAA CATCGGGGCG GA 
(2) INFORMATION FOR SEQ ID NO: 7: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 33 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER" 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 7: 

ATCGGGGCGG ACGGTTGTTG GGTGTCGGGC GCG 

(2) INFORMATION FOR SEQ ID NO: 8: 

<i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 25 base pairs 
<B) TYPE: nucleic acid 
<C) STRANDEDNESS: single 
(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc - "PRIMER" 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:8: 
GAGTATCGTG TGTACTGGCG GCACC 
(2) INFORMATION FOR SEQ ID NO: 9: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER" 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 9: 
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GGTCTCGGTG AGCCCAGGTT CAATGTCGAT 
(2) INFORMATION FOR SEQ ID NO: 10: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER* 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 10: 
GGTCTCGGTG AGCCCATGTT CAATGTCGAT 
(2) INFORMATION FOR SEQ ID NO: 11: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 27 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
<D> TOPOLOGY : linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER" 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 11: 
GAGGACACGT ACTGGAACGG * CAACCCG 
(2) INFORMATION FOR SEQ ID NO: 12: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 60 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 
(A) DESCRIPTION: /desc = "PRIMER" 



(xi) SEQUENCE DESCRIPTION : SEQ ID NO: 12: 
TACCCTGAGG ACACGTACAA CGGCAACGGC AACTCGCAGG GCAACCCGTG 
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CLAIMS 

What is claimed is: 

1. A fungal glucoamylase including a mutation 
pair Asn20Cys coupled with Ala27Cys forming a disulfide 
bond between the two members of the pair. 

2 . The glucoamylase as set forth in claim 1 
wherein the mutation provides increased thermal 
stability and reduced isomaltose formation. 

3 . The fungal glucoamylase as set forth in claim 
1 and including at least one mutations selected from 
Table 13 wherein cumulative thermal stability is 
provided by the additional mutations. 

4 . The fungal glucoamylase as set forth in claim 
1 further including mutations Ser30Pro, Glyl37Ala 
wherein cumulative thermal stability is provided by the 
additional mutations. 

5 . The fungal glucoamylase as set forth in claim 
1 and including at least one mutation from Table 14 
wherein cumulative reduced isomaltose formation is 
provided by the additional mutations. 

6. The fungal glucoamylase as set forth in claim 
1 further including a 311-314Loop mutation wherein 
cumulative reduced isomaltose formation is provided by 
the mutation. 

7. A fungal glucoamylase including a 311-314Loop 
mutation. 

8. The glucoamylase as set forth in claim 7 
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wherein reduced isomaltose formation is provided by the 
mutation, 

9. The fungal glucoamylase as set forth in claim 
7 and including at least one mutation from Table 14 
wherein cumulative reduced isomaltose formation is 
provided by the additional mutation. 

10 . A fungal glucoamylase including a mutation 
Ser411Ala. 

11 . The glucoamylase as set forth in claim 10 
wherein increased pH optimum and reduced isomaltose 
formation is provided by the mutation . 

12 . The fungal glucoamylase as set forth in claim 
10 and including at least one mutation from Table 15 
wherein cumulative increased pH optimum is provided by 
the mutations. 

13 . The fungal glucoamylase as set forth in claim 
10 and including at least one mutation from Table 14 
wherein cumulative reduced isomaltose formation is 
provided by the mutations. 

14. A fungal glucoamylase including a Ser411Ala 
mutation and a mutation pair Asn20Cys coupled with 
Ala27Cys forming a disulfide bond between the two 
members of the pair. 

15. The glucoamylase as set forth in claim 14 
wherein increased thermal stability, increased pH 
optimum and reduced isomaltose formation are provided 
by the mutations. 

16 . A fungal glucoamylase including a Ser411Ala 
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mutation and a mutation pair Asn20Cys coupled with 
Ala27Cys forming a disulfide bond between the two 
members of the pair and a 311-314Loop mutation. 

17. The glucoamylase as set forth in claim 16 
wherein increased thermal stability, increased pH 
optimum and reduced isomaltose formation are provided 
by the mutations. 

18 . A vector containing the cDNA for an 
engineered glucoamylase as set forth in claims 1-17. 

19. A host cell transformed with the vector of 
claim 18. 

20. A fungal glucoamylase as set forth in claims 
1-17 wherein the glucoamylase is an Aspergillus 
glucoamylase . 

21. The glucoamylase as set forth in claim 20 
wherein the glucoamylase is Aspergillus awamori 
glucoamylase. 

22. A method to obtain a fungal glucoamylase with 
decreased thermal inactivation by designing mutations 
having decreased conformational entropy of unfolding, 
or increased stability of a-helices, or increased 
disulfide bonds or hydrogen bonding and electrostatic 
interactions and hydrophic interactions and Vanderwalls 
interactions and packing compactness. 

23 . A method to obtain a fungal glucoamylase with 
reduced isomaltose formation by designing mutations to 
decrease the a- (1, 6) -glucosidic linkage affinity. 

24. A method to obtain a fungal glucoamylase with 
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increased pH optimum by changing the polarity, charge 
distribution and hydrogen bonding in the 
microenvironment of the catalytic base Glu400. 

25. A method of selecting mutations for fungal 
glucoamylase to be used for constructing glucoamylases 
with cumulative mutations by 

designing and generating individual mutations by 
site directed mutagenesis; 

screening the individual mutations and selecting 
those which show at least increased pH optimum, 
decreased irreversible thermal inactivation rates or 
reduced isomaltose formation; 

performing site directed mutagenesis to produce 
enzymes carrying at least two of the isolated selected 
mutations for either increased pH optimum, decreased 
irreversible thermal inactivation rates or reduced 
isomaltose formation; and 

screening for cumulatively additive effects of the 
mutations on pH optimum, thermal stabilizing or reduced 
isomaltose formation by the produced enzymes carrying 
at least two of the isolated selected mutations. 
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